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Abstract 

Motivated by the study of non-abelian Chern Simons vortices of non-topological type 
in Gauge Field Theory, see e.g. [33J 33),[21], we analyse the solvability of the following 
(normalised) Liouville-type system in presence of singular sources: 

! — Aui = e Ul — re“ 2 — ANttSq, 

—A 1 L 2 = e“ 2 — re Ul , 

= h eUl and ft = sr /« 2 e “ 2 - 

with r > 0 and N > 0. 

We identify necessary and sufficient conditions on the parameter r and the ’’flux” pair: 

(Pi, ^ 2 ), which ensure the radial solvability of (1) T . 

Since for r = p problem (1) T reduces to the (integrable) 2X2 Toda system, in particular 
we recover the existence result of m and m, concerning this case. 

Our method relies on a blow-up analysis for solutions of (1) T , which (even in the radial 
setting) takes new turns compared to the single equation case. 

We mention that our approach permits to handle also the non-symmetric case, where in 
each of the two equations in (1) T , the parameter r is replaced by two different parameters 
7~i > 0 and T 2 > 0 respectively, and when also the second equation in (1) T includes a 
Dirac measure supported at the origin. 
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1 Introduction 


In recent years there has been a growing interest towards the understanding of non-abelian 
vortex configurations supported by Chern-Simons Gauge Field Theories, also in view of their 
connection to the delicate issue of ’’monopole confinement”. See [331 ,3lj and references 
therein. 

In this context, a successful way to detect vortices is to identify the BPS-sector of the theory, 
since in such regime vortex configurations simply correspond to (static) solutions of the 
so called self-dual equations of Bogomolnyi type, and saturate the minimal energy allowed 
by the system. For this reason, it has been useful to invoke ’’duality” and formulate the 
theory within the general framework of J\f = 2 Supersymmetric (SUSY) Field Theory, in this 
direction see for example: [27U281 EHES (Ml EH EO] for more details . 

We observe that, when the Chern-Simons Lagrangian is taken into account then the theory 
can attain self-duality only with the help of a six-order scalar potential field, see [261 Ell El], 
in place of the more familiar quadratic (double-well) potential of the Maxwell-Higgs model, 
see [30] • As a consequence, the system acquires additional vacua states, in the sense that 
now, the unbroken vacuum must coexist together with broken and partially broken vacua 
states. In turn, this feature give rise to new classes of self-dual vortices, each supported 
(asymptotically) by different vacua states. 

On the other hand, in the non-abelian setting, it is quite difficult to handle self-dual vortices 
with mathematical rigour, since even the reduced first order self-dual equations present serious 
analytical difficulties yet to be overcome. Nonetheless, it is possible to focus on special classes 
of vortices satisfying some useful and physically consistent ansatz, which permit to formulate 
the self-dual Chern-Simons vortex equations into a form suitable for the approach in [30] to 
apply. In this way, one is lead to the so called ’’Master equations”, which govern such class 
of vortex configurations, and take the form of planar elliptic systems involving exponential 
nonlinearities and involve Dirac measures supported at the vortex points, see for example 

OS 113133[IISES ED [26]. 

By means of the Master equations and their variational structure, it is possible to construct 
planar topological vortices (asymptotically gauge equivalent to the unbroken vacuum) by a 
minimisation procedure, see [THESES]. Similarly, after the variational approach introduced 
in [SEE! for the abelian case, it is possible to use variational methods also to establish multi¬ 
ple periodic non-abelian Chern-Simons vortices for various theories, see e.g. [MIESESES. 

On the contrary, the study of planar non-topological vortices, asymptotically gauge equivalent 
to the broken vacuum, still remains unsatisfactory, in spite of the rather complete description 
available for the abelian case, see |60| eh is Eg. Indeed, even when all the vortex points 
are superimposed, say at the origin, and one seeks radially symmetric vortices (about the 
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origin), the corresponding radial differential vortex problem turns out to be quite involved. 
In this direction, we mention the contributions in m, Esi, m and }18j which provide useful 
general information about the nature of entire radial vortex solutions, and reveals their rich 
structure compared with the abelian situation. 

On the basis of such results and as a way to progress further towards the understanding of 
non-topological vortices, here we propose to adopt the ’’perturbation” strategy introduced 
by Chae-Imanuvilov in [8] for the abelian case, see also mm for non-abelian settings. 
Actually, to carry out such ’’perturbation” approach, the first important step is to provide 
rather accurate informations about the solvability of Liouville - systems involving Dirac mea¬ 
sures, which occur as ’’‘limiting’” problems in the perturbation argument. 

Interestingly, the relevance of such class of Liouville systems has emerged already in several 
other contexts, see [131031051061 HU ED S3] and references therein. Their study has concerned 
mainly the so called ’’cooperative” case, where all the entries of the coupling matrix are 
assumed positive, and we refer to [121 13, ED E3, El ED, E2J S3) 1331 EH, ED], for rather 
complete results about this situation. 

However, for the non-abelian models considered here, the corresponding Liouville systems are 
of ” competitive” - type, and their solvability is far from understood, except for the very special 
case of the integrable Toda system. We recall that, the (m x m)-Toda system is characterised 
by its coupling matrix being the Cartan matrix of the (gauge) group SU{m + 1), and it 
describes vortices for a relativistic Chern-Simons model, see [26]. 

Lin-Wei-Ye in [50] have identified all the (conformal) invariances enjoyed by the ’’singular” 
Toda system, and used them to obtain a full description of the corresponding solution set and 
of their linearised problem. The classification result in [50] was motivated and inspired by 
the analogous classification result obtained by Jost-Wang [43] for the ’’regular” Toda system, 
where the Dirac measures were neglected. 

Also note that for m = 1 (corresponding to the abelian case discussed in [8]), the result in 
[50 ] just reduces to the well known classification result for solutions of ’’singular” Liouville 
equations obtained in m (see also PEE]), and corresponding linearised problem discussed 
in [21]. 

The complete analysis provided by in |5D] for the singular Toda system has allowed Ao-Lin- 
Wei in mm to pursue the ’’‘perturbation” approach of [8], and construct non-topological 
vortices for the Chern-Simons model proposed by Dunne in [ 26] . with any gauge group of 
rank 2. 

The main purpose of this paper is to analyse the solvability of (not necessarily integrable) 
planar ’’‘singular’” Liouville systems, as they occur in various Chern-Simons vortex problems, 
see e.g. [Ml SUES], or in other physical context decribed for example in @11051 06] 03, and 
which include the Toda system as a particular case. 
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In our situation, the systems enjoys only a scale invariance property, and so the admissible 
” flux-pair” corresponding to a solution is no longer a fixed value (as in the conformal / 
integrable case) but now it describes a curve, which we will be able to describe completely in 
the radially symmetric case. 

More precisely, we provide necessary and sufficient conditions for the radial solvability of a 
general singular 2X2 Liouville system of Toda type, in terms of the ’’flux pair” of a solution 
(given by certain integral values), which relates to the total energy carried by the correspond¬ 
ing vortex, see Theorem 14. 121 for the precise statement. In certain cases, our results are sharp 
also when the radial assumption is dropped, see Corollary 14.141 . 

In certain sense, our analysis should be considered the continuation of that in (55] (see also 
[55] ) of’’cooperative” systems. 

To simplify notation and to better clarify the arguments involved, we shall work mainly in 
the symmetric case, although it will be clear by our discussion how to modify the given 
arguments in order to treat non-symmetric systems as well. We mention also that, [29] , [20] 
and [32] contain some particular cases of our result. 

A central aspect of our proof is to provide an accurate blow up analysis for solutions of 
problem (1) T , which relies in a crucial way upon the radial symmetric assumption. In fact, to 
obtain a similar blow-up description for general (non-radial) solutions would be much harder 
to attain, as we can see already for the integrable Toda system analysed first in [12] and more 
recently in m (see also references therein), or for a degenerate system given by a Cosmic 
Strings equation discussed in [55] , 

In the next section, we provide the necessary preliminaries and the precise statements of our 
results. 

2 Preliminaries and Statement of the Main Results 

Motivated by the delicate issue of ’’‘quark confinement”’, Gudnason in [33j 1H33 introduced 
a non-abelian Chern-Simons model formulated within a JV = 2 Supersymmetric (SUSY) 
Field Theory, with a general gauge group of the type: G = U(l) x G' allowing solutions 
with orientational modes. For the model in M EH, the author identifies the BPS-sector 
of the theory and the corresponding self-dual equations. In particular, when G’ = 50(2) 
or G' = US p ( 2), Gudnason in [ [33! 133] introduced some meaningful physical ansatz on the 
structure of the vortex solutions, by which (as in [40]) the corresponding self-dual equations 
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reduced to the following set of Master’s equations: 

J A u = a 2 (e u+v + e u ~ v - £)(e u+v + e u ~ v ) + aP(e u+v - e u ~ v ) 2 + S qi j 

[Av = a/3(e u+v + e u ~ v - £)(e u+,; - e u ~ v ) + p 2 (e^ u+v ^ - e 2 ^)) 2 + 4ttY]Ii <W 

( 2 . 1 ) 

with a = and f3 = ^ and where fci and k 2 are respectively the coupling constants of the 
17(1) and the G' part of the gauge field; moreover £ > 0 is known as the Fayet-Iliopoulus 
parameter. 

The set of points (repeated with multiplicity): 

Si = Ui,j,j = l,-,rcj} 7 = 1,2 (2.2) 

where the Dirac measures are supported, are known as the vortex points, and they carry a 
relevant physical meaning, see [33l f3T j. In order to avoid singularities in the component u — v, 
we need to require that, 


ri 2 < ri\ and S 2 C S\. (2-3) 

For more details about the model and the derivation of (12.11) we refer to [33) SHJ and [35]. 
Next, we observe that we can express (El]), in a form more familiar within the context of 
self-dual non-abelian Chern-Simons vortices, see 123 ESI ED. 

To this purpose, we let, 


A - 7—) 2 and = ^ = h 
V 2) a k 2 

and we introduce the new unknowns: 


(2.4) 


U\ = U + V , U2 = U ~ V 

in terms of which, the system ED takes the form: 

Arq + A 2 (y!= 2 k i,i eUi ( x - Ej= 2 hi^ 3 )) = 4?r E?=i <W 
Ari 2 + A 2 (EL2 k ^ Ui ( x - Ej = 2 hj eUj )) = 47r E^i 
with coupling matrix K = (kij) given as follows: 

1 / 1 + k 1 -k 
~ 2 ^ 1 - k 1 + k 

and (by recalling (12.21) 1: 

Ni = m + n 2 and {p 1>j ,j = l,..,N 1 }:=S 1 llS 2 , 


(2.5) 


( 2 . 6 ) 
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N -2 = ri\ — n 2 and {p 2 ,j,j = 1, -,N 2 } := <51 \ S 2 . 

Observe that, for k = 3, the coupling matrix (12.61) coincides exactly with the Cartan matrix 
of the group SU(3) given as follows: 


K = 


2 

-1 



(2.7) 


Therefore, in this case, (12.51) reduces to the well known (integrable) 2 X 2 - Toda system, 
whose solutions describes .SC/^-vortices for the relativistic Chern-Simons- model proposed 
by Dunne [26). More generally, the model in [26] is formulated in terms of any general (semi¬ 
simple) gauge group G, with A = and /1 > 0 the Chern-Simons coupling parameter, 

see also [ZD EH] for additional details. 

In the more general framework of the J\T = 2 — SUSY, we mention also the model proposed in 
[55] concerning a non-abelian Chern-Simons theory with ’’flavors” and gauge group: 17(1) x 
SU(N), where a system of the type (12.51) occurs as the Master equations governing the 
corresponding self-dual vortex configurations with coupling matrix: 


l / N — 1 + k 1 - k 

N y (N — 1)(1 — k) 1 + (N — l)k 


( 2 . 8 ) 


where again k = and k\ > 0, k 2 > 0 are respectively the coupling constants of the 17(1) 
and the SU(N) part of the gauge fields, and A = ( 9 ^) • 

Notice that, in all the given examples, and actually in other examples of physical interest, 
the corresponding coupling matrix K is characterised by the following property: 


k lyi > 0 and k ltl + k± :2 = /c 2 ,i + k 2 ,2 = 1. (2.9) 

To construct planar vortices, we seek solutions of (12.51) in K 2 , satisfying the following inte- 
grability condition: 

ki,ie Ui 1^1 - kije Uj 

which ensures that, the corresponding vortex admits finite energy , see [m [331 1E1 EH]- 
By keeping in mind (12.91) . from (12.101) we identify the following suitable boundary conditions, 
as Ixl —y 00 : 


X 


^ € T 1 (M 2 ), 1 = 1,2; (2.10) 


topological: 

Ui(x) —> 0 i = 1,2, 

(2.11) 

non-topological: 

Ui(x) —>• —00 i = 1,2, 

(2.12) 
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mixed type: 


(2.13) 


mi(.t) —>• —oo and U 2 (x) —>• — log ^ 2 , 2 ? 
or 

U‘ 2 (x) —>• —00 and u±(x) —>• — logfc^i. 

From the physical point of view, we mention that, topological solutions give rise to vortices 
asymptotically (at infinity) gauge equivalent to the unbroken vacuum state, non-topological 
solutions give rise to vortices asymptotically gauge equivalent to the broken vacuum state, 
and finally mixed-type solutions allow for vortices asymptotically gauge equivalent to par¬ 
tially broken vacua, see mi ED ED for details. 

To further substantiate the physical pertinence of the boundary conditions specified above, we 
observe that (when (12.91) holds) every solution of (12.51) satisfying either one of the boundary 
conditions (I2.11|) . (I2.12|) and (I2.13|) . verifies: itj < 0 in M 2 for every i = 1, 2, consistently with 
the physical applications. 

Furthermore, in the case of interest here, namely: 

ki t i > 0 and h h j < 0 with z / j = 1,2, 

and when p t j = 0 for every l = 1,..., iVj and z = 1,2, then Huang-Lin in (3?, 38] have shown 
that every entire radially symmetric solution of (12.51) must satisfy either one of the boundary 
conditions above, and in particular the integrability condition (12.101) automatically holds in 
this case. 


The construction of topological solutions for (12.51) has been carried out first by Yang in ED, 
and more recently in [55] and m for more general choices of coupling matrices including 
(12.61) and (12.81) as particular cases. We also mention that the existence of multiple periodic 
solutions for (12.51) was established first by Nolasco-Tarantello in [56] for the choice of coupling 
matrix (ED), and more recently in [36], [35] and m for more general choices including (12.61) 
and (|2.8I) respectively. 

In this paper, we focus on the construction of non-topological solutions for (12.51) . 

To this purpose, we recall what happens for the embedded abelian configurations, corre¬ 
sponding to the situation where, 


and 


N = N\ = N 2 and (after relabelling) p\j = p 2 j = pj, j = 1,.., N, 


U\ = U 2 = U 


with u a solution of the abelian Chern-Simons vortex equation (cfr. [39l 0TJ) : 
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N 

Au + Xe u (l — e u ) = 47r ^ 5 Pj in M 2 , (2.14) 

l=i 

(recall (12.9p ). For a detailed study of (12.141) we refer to [68]. 

In particular, we know that a non-topological abelian vortex will correspond to a solution u 
of (|2.14p satisfying: e u € L 1 (M 2 ). To construct such class of solutions, we let: 

u £ (x) = u + 2 log i, e > 0; 

and observe that u £ satisfies: 


N 

— A u £ = Xe Ue — Xe 2 e 2us — 47r ^ d £Pj in M 2 , 

l=i 


(2.15) 


and 



(2.16) 


Therefore, (at least formally) non-topological solutions of (12.141) can be sought as ’’pertur¬ 
bation” of solutions of the following ’’limiting” problem, obtained by letting e —>• 0 in (12.151) . 
namely: 

- Au = Xe u - 4 itN5 0 in M 2 , (2.17) 


and, 

e u £ L 1 (M 2 ). (2.18) 

Solutions of (12.171) . (12.181) are explicitly known (see ([2.291) below, and [57], [HJ HZ]), and 
they were used by Chae-Imanuvilov in [8] together with suitable ’’perturbation” techniques 
to obtain a family of non-topological solutions for (I2.14p . we refer to [8] and [67] fo details. 
Actually, the subsequent analysis in m and [22] has provided a rather complete description 
about the whole set of non-topological solutions for the single equation (12.141) . 


Unfortunately, the analysis of non-topological solutions for the non-abelian vortex system 
(12.51) is far from complete, in spite of the results in [1, [2], 3?) 38l 123] . 

For example, one we can perform a similar scaling for a solution (u±,U 2 ) of (12.51) as follows: 

ui i£ (.x) = ui + 2 log i , u 2 , £ (x) = u 2 + 2 log i (2.19) 

and realise that, in order to pursue a ’’‘perturbation’” approach similar to [8], it is important 
to understand the complete solvability of the following singular Liouville system: 
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( 2 . 20 ) 


— Au\ = ki t ie Ul + /ci )2 e“ 2 — 4nNi5o 
< -Au 2 = k2,ie Ul + k 2) 2e U2 — 4'kN 2 5o 
e ul ,e U2 € L 1 (R 2 ) 

which ’’formally” we derive by letting e —> 0 in the system satisfied by (ui j£ ,U 2 j£ )- 

The analysis of (12.20p is interesting in its own, since similar Liouville systems (in absence 
of Dirac sources, i.e. N\ = N 2 = 0) have emerged in other contexts, such for example 
in statistical mechanics, see EH S3 06) HZ] and nana, population dynamics and plasma 
physics, see eh m Ea and references therein. 

More precisely, concerning (12.201) . we address the question of identifying the pairs: (/3 i,/ 3 2 ) 
for which problem (12.201) admits a solution satisfying: 

Pi = -}- [ e Ul dx, fh = ^~ [ e“ 2 dx. (2.21) 

2tt Jr 2 27T Jr 2 

This is a natural question to investigate also from the physical point of view, since the values 
[3 1 and /3 2 indicate the energy level of the corresponding vortex configuration. 

We also notice that, in the context of radial solutions, problem (I2.20|) . (12.211) can be turned 
conveniently into a more standard boundary value problem. 

Actually, in the radial setting, we have a complete answer about the question of solvability 
for the so called cooperative systems, even for general m x m — systems , where the coupling 
matrix is assumed to satisfy the following: 

K = (ki j)i j = i m is symmetric, irreducible with non-negative entries. (2.22) 

Solutions of cooperative systems in absence of Dirac measures are always radially symmetric, 
see L2U •, and they have been investigated in [12) OGQ 2L 52) 153) 63) [6D [70, 12) 03], also in 
connection with sharp versions for systems of the Moser-Trudinger inequality, see m , m , 
or via duality, by invoking sharp versions of the log — HLS inequality for the Free Energy, 

(see eh m E]). 

More recently, those results have been extended by Poliakovsky-Tarantello in [58, [50], to 
describe the sharp radial solvability for ’’singular” Liouville systems including the Dirac mea¬ 
sures and with coupling matrix K satisfying (12.221) . 

It should be noticed that in this case, non-radial solutions do occur, as exhibited in [58]. Fur¬ 
thermore, [58] EH] cover also the degenerate case, (unlike El] an d [52)153]) which includes a 
class of Liouville-type equations arising (as a degenerate 2x2-system) in the study of self-dual 


9 


















Cosmic Strings, see [58], [59] and [69] for details. 


For problem (12.201) . (12.211) the results [58] [59] can be summarised as follows: 

Theorem 2.1 ([T2, 151 [211 59]). Assume that, 

h,j = kj,i + 0, for ijtje {1, 2}, (2.23) 

and let Ni > — 1, i = 1,2. Then the following conditions: 

h,ifl + k 2 , 2 Pl + 2fci, 2 0i0 2 - 4(JVi + l)ft - 4(!V 2 + 1)0 2 = 0 (2.24) 

+ <0, * = 1,2 (2.25) 

are necessary for the existence of a radial solution to (|2.20l) . (12.211) . 

Furthermore, if (Cooperative Systems) 

h,i > 0 , kij = k hi >0, i / j € {1,2} (2.26) 

then (12.241) . (12.251) are also sufficient for the radial solvability of (12.201) . (12.211) . 

In addition, if —1 < N < 0 then every solution of (|2.20[) . (12.211) is indeed radially symmetric. 

□ 

Remark 2.2. Observe that, the case &q )2 = k 2i \ = 0 is not interesting, as (12.20(1 decouples 
into two singular Liouville equations of the type (12.171) . and it is solvable if and only if 
ki,il3i = 4(lVj + 1), i = 1,2. 

The conditions: (I2.24|) - (I2.25|) rely in a crucial way on the scale invariance of problem (12.201) . 
(12.211) under the following transformation: 

u it \(x) = Ui(Xx) + 2 log A, * = 1,2 and A > 0; (2.27) 

in the sense that, 

(ui,u 2 ) satisfies (12.201) . (12.211) (ui t \,U 2 ,\) satisfies (12.20(1 . (12.211) : V A > 0. 

The following uniqueness result holds: 

Theorem 2.3 (|52] [59]). Assume ([2.26D . then for every (0i,0 2 ) satisfying (|2.24j) . (12.251) 
there correspond a unique radial solution of (12.201) . (12.211) modulo the transformation (12.271) . 
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The above uniqueness result is based upon a non-degeneracy property for the linearised prob¬ 
lem around a given radial solution, in the sense that (in the radial framework) it admits as 
only degeneracies those originated by the scale invariance (12.2711 : we refer to [52 4 [59] for 
details. This fact extends what we know to happen for the decoupled system (12.2011 . with 
fci,2 = k 2 ,i = 0. 

Clearly, Theorem 12.11 and Theorem 12.31 apply to Gudnason’ s model (with coupling matrix 
(12.611 and gauge group G = U(l) x G'), when the strength k 2 of the G"-component of the 
gauge field overpowers the strength k\ of the 17(l)-component, that is: k = < 1. So, on 

the basis of the uniqueness (and non-degeneracy) property of radial solutions given by The¬ 
orem 12.31 in this case it should be possible to use bifurcation techniques in order to obtain 
non-topological (radial) solutions for the whole system (12.511 . 


As already observed, in case N\ > 0 or N 2 > 0, then we no longer expect every solution 
of (12. 201) . (12.2111 to be radially symmetric. Such a symmetry breaking phenomenon already 
occurs for the single Liouville equation (12.1711 - (12.181) . where every solution u must satisfy: 


/ 

Jr- 


e u dx = 4(N + 1), 
and (in complex notation) takes the following form: 


(2.28) 


u(z) = log 


8{N + 1) 2 h\z 


2N 


(2.29) 


_A(l + /r|z Ar + 1 + 6| 2 ) 2 _ 

with [i > 0 and b £ C satisfying: b ^ 0 if and only if N € N (see [57]). In particular, we see 
that only for N £ N, all solutions to (12.1711 - (12.181) are radially symmetric. 


On the other hand, it seems reasonable to expect that (12.241) . (12.251) still provide sharp nec¬ 
essary and sufficient conditions for the solvability of (12.201) . (12.211) . also beyond the radial 
situation. However, while the Pohozaev-type identity ([2.2411 remains valid for non-radial solu¬ 
tions of (I2.20I) - (I2.21D . as it relies only upon the scaling invariance (12.271) and the asymptotic 
behaviour of solutions at infinity (see Proposition 12.41 below), it is less obvious to verify 
whether the condition (12.2511 holds also for non-radial solutions. 

More precisely, by letting: 


the following holds: 


— hi,i/3i + ki,jPj, if - J € {!> 2} 


(2.30) 


Proposition 2.4. Let Ui € L\ oc { R 2 ), i = 1,2, and suppose that (u\,u 2 ) satisfies (I2.2U|) . 
(12.2111 (in the sense of distributions) for a given pair (/3 i,/ 3 2 ). We have: 
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(2.31) 

(2.32) 


(i) 


Ui(x) = -( Pf° + 2Ni) log |x| + 0(1), for \x\ > 1 

rd r Ui(x) — »■ — (/3f° + 2 Ni), deUi(x) — >• 0 as r = |x| — >• +oo 

with /3f° given in (|2.3U|) . * = 1,2; and (r,9) the polar coordinates in M 2 . 

In particular, there holds: 

P? = ki,iPi + hjPj > 2 (Ni + 1), Vi + j € {1,2}. (2.33) 

(ii) If K = 1,2 satisfies: 

hi ,2 • ^ 2,1 > 0, (2.34) 

then the pair (/3 i,@ 2) verifies the following identity: 

h,i\k 2 ,i\Pl + k 2:2 \k 1;2 \Pi + 2fc lj2 1 k 2 ^i\(3 1 p 2 - 4(N X + 1 )\k 2 ,i\Pi - 4(JV 2 + l)\k h2 \p 2 = 0. (2.35) 

□ 

Clearly (12.351) reduces to (I2.24|) when K is symmetric, and in fact in such case, Proposition 
IP is essentially established in mmm- For completeness we shall indicate its proof in 
the following section. 

Furthermore, if one of the off diagonal entries of K vanishes, (so that (12.341) holds with equal 
sign), then (12.35[) simply gives back what we already know about the single Liouville equation 
(12.171) . Actually in this case, one is left to analyse a single planar equation, which can be 
casted as a weighted Liouville equation and in this way it can be handled as in [49], [58] . [25]. 

So the true interesting case to investigate occurs when k\ i2 ■ k 2j \ > 0. 


Remark 2.5. Condition (12.331) does not appear for the cooperative system (12.261) . since in 
this case it is a consequence of (12.241) and (12.251) . Indeed if by contradiction we assume for 
example that, 


ki^iPi + Aq, 2 P 2 < 2(lVi + 1), 

then 2k\,i[3'f + 2k\. 2 fi\fi 2 < 4(lVi + 1)P\ and from (12.241) we see that necessarily k 2)2 p 2 > 
4( N 2 + 1), in contradiction with (12.251) . 

The main focus of this paper is to investigate the solvability of (I2.20l) - (|2.2ip in the more 
delicate ’’competitive” situation, where we assume: 

k{ t i > 0 and kij <0 for i j £ {1,2} (2.36) 


12 































as it occurs for (12.61) when k = > 1. 

In particular (12.361) covers the 2X2 -Toda system (12.71) . which is the only situation where, 
so far, we can describe the whole solution set of (12.201) . Indeed, it has been shown in [50] . 
that in this case the system gains full conformal invariance and, in analogy with the single 
Liouville equation (12.171) . (cfr. |5T]), it is possible to identify explicitly the corresponding 
solutions. In particular, we may conclude the following about the values (5\ and /?2 in (I2.2ip : 


(12.201) — (12. 2 II) with K in (12.71) is solvable if and only if f5\ = /3 2 = 2(A r i + 1V 2 + 2). (2.37) 

We refer to [03| and [50] for a detailed description of the Toda system and its relevant con¬ 
served quantities. 

Under the assumption (12.361) . the necessary condition ([2.331) plays a central role towards the 
solvability of (I2.20I) - (I2.21I) and it implies the following: 

Corollary 2.6. If ki } j < 0 for i ^ j £ {1,2}, then 

ki t i > 0, is {1,2} and det(K)> 0, (2.38) 

are necessary conditions for the solvability of (12.201) , (12.211) . In particular in the symmetric 
case, the solvability of (12.201) , (12.211) requires that K is strictly positive definite. 

From (12.381) . we see in particular that (12.361) and (12.261) are truly complementary and fully 
account for (12.231) . 

In order to proceed further, we observe that under the assumption ()2.35j) we can make an 
useful normalisation, simply by setting: 

Ui(x ) = Vi(x ) + 2 Ni log |x| - log k it i, i = 1, 2; (2.39) 

and in terms of the new unknowns ( , ri,'r 2 ), we obtain the following problem: 

—Avi = \x\ 2Nl e Vl — ti\x\ 2N2 e V2 

< —Av 2 = \x\ 2N2 e V2 — T 2 \x\ 2Nl e vi (2.40) 

h Jr2 \x\ 2Ni e Vi = /3ik it i, * = 1,2 

with, 

ti = > 0 and t 2 = — > 0. (2-41) 

« 2,2 «11 
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Remark 2.7. For the coupling matrix given in m, problem (12.401) occurs with, 


n=r 2 = > °> f° r k > L ( 2 - 42 ) 

In particular, the singular Toda system (where k = 3) leads to (12.401) with, 

ti = t 2 = \ and [ \x\ 2Nl e Vl = [ \x\ 2N2 e V2 = 4(JVi + N 2 + 2), (2.43) 

2 2vr J K 2 2tt J r2 

see (1071) and (12.401) . 

To simplify notations and to minimise technicalities, in this paper we shall focus on the study 
of ()2.401) in the symmetric case, where we assume: 

Ti = t 2 := t > 0, (2.44) 

and we shall give indications on how to deal with the non-symmetric case: ti ^ t 2 and 
t\t 2 > 0, which will be discussed in details in a forthcoming paper. 

For the original system (I2.20p , the validity of (12.441) amounts to assume the following stronger 
version of (|2.35l) : 

ki.i > 0 and for i/j S {1,2}. (2-45) 

We mention that, a system of the type (12.201) for which (12.451) holds, is referred to as a 
collaborating system, see [63, 611. 

Since our original motivation was to analyse Gudnason’ s model, w.l.o.g. we can further 
assume N\ > N 2 \ and again to reduce useless technicalities, we shall specify the following: 


N\ = N > 0 and N 2 = 0 (2-46) 

Therefore, from now on we shall focus on the solvability of the following problem: 




—Av\ = \x 


_ \ry\2N Vl 


e 1 — re 


V2 


2N p vi 


—Av 2 = e V2 — t\x 


Pi = h /r2 \x\ 2 N e n , P2 = ^ f R 2 


with r > 0 and N > 0. 

According to the discussion above, already we may claim the following about the solvability 
of (P) T : 


Corollary 2.8. If t > 0 and > 0, i = 1,2, then the following conditions are necessary for 
the solvability of (P) T : 
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re ( 0 , 1 ) 

Pi + Pi - 2 T p x p 2 - 4(N + 1 )p! - 4p 2 = 0, 

(2.47) 

(2.48) 

Pi — t/3 2 > 2 (N + 1) and 

P2 - rPi > 2 . 

(2.49) 

In addition the condition: 

Pi > 4(N + 1) and 

P 2 >4, 

(2.50) 

is necessary for the existence of a radial solution. 


□ 


In contrast to the case where r < 0 (see part b) of Theorem 12.ip . we do not expect in general 
that conditions (|2.47|) - (|2.49j) (and (12.501) ) are also sufficient for the (radial) solvability of (P) T . 
Indeed, it suffices to recall what happens for the Toda system, where by (12.431) and (12.461) we 
have: 

if t = i then Pi = p 2 = 4(IV + 2), (2.51) 

and we see clearly that (12.511) is a much more stringent condition than (I2.47D - (I2.50I) . 

On the other hand, by virtue of part a) in Theorem 12.11 which applies to (P) r , and by keeping 
in mind part b) (which applies to (P) T only when r < 0), by way of continuity, we expect 
that (12.481) and (12.501) should still provide necessary and sufficient conditions for the radial 
solvability of (P) T , when r > 0 is sufficiently small. 

We establish this fact and prove the following: 

Theorem 2.9. Let N > 0, then for any, 


0 < r < 


:= r 0 (JV) 


JV + l + ^(fV+l) 2 +4 
conditions (12.481) and (12.501) are necessary and sufficient for the radial solvability of (P) 


(2.52) 


□ 

Clearly to (IV) < \ (recall that r = \ corresponds to the Toda system) and, as we shall see, it 
corresponds to the sharp value for which conditions (|2.48j) . ( 12 . 501 ) imply (|2.49l) . consistently 
with Remark 12.51 

By combining Theorem 12.II and 12.91 we can conclude the following about problem (P) T . 
Corollary 2.10. For every 0 7 - r < tq(N) the conditions (12.481) and 


t ^ t (/3i-4(JV + 1))<0 ; ^(/3 2 -4)<0 (2.53) 

|r| \t\ 

are necessary and sufficient for the radial solvability of (P) r . 
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□ 

We suspect that, also the uniqueness (and non-degeneracy) result of Theorem 12.31 should 
remain valid for radial solutions of (P) T , whenever r < to (IV), and not only when r < 0 . 

Interestingly, for r sufficiently ’’close” to 1, the conditions (12.481) and (12.491) are the ones 
responsible for the solvability of ( P) T . 

To be more precise we consider the following polynomial, whose role will become clear in the 
following, 


fii (r) = 2(1 - 2t 2 )(1 + 2 t(N + !))-! 


and by direct inspection we check that: 


(r) admits a unique positive zero, denoted by t\ = t\(N), 

a,,dTie Q-7l) 


(2.54) 


Theorem 2.11. For any t € [ti, 1 ), with n = t\{N ) defined in (|2.54p . we have that the 
conditions (12.481) and (12.491) are necessary and sufficient for the solvability of (P) T . 


□ 

Actually, under the assumptions of Theorem 12.111 we shall see that the conditions (12.481) and 
(12.491) imply (12.501) and permit to establish the existence of a radial solution for ( P ) T . Hence, 
it is tempting to conjecture that ( P) T may admit only radial solutions in this case. 


The case where r € ( tq,t± ) is more delicate, as it includes the Toda system (r = 1/2), where 
we know that the necessary conditions (12.481) - (12.501) are no longer sufficient for solvability, 
see (12.511) . 

In fact, for t € (to,ti) we are able to identify new necessary and sufficient conditions for the 
(radial) solvability of ( P ) T , which indeed reduce to (12.511) when r = 1/2. 

To this end we observe that the conditions (I2.48P and (12.491) define a graph in the first 
quadrant of the (ffi, / 32 )-plane. 

More precisely, for r € (0,1), we let: 


(r) = —^ (N + 1 + r + t \J(N + l) 2 + 2t(N + 1) + l') 

2 T / x 7 (2 ’ 55) 

/?i := /3 1 (r) = ^—2 (N + 1+T+ ^{N+1) 2 + 2t(N + 1) + 1) > 4 {N + 1) 

so that, /3 1 (r) < and similarly, we let: 
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(2.56) 


P_ 2 ■■= P 2 (.t) = YZ^2 ( X + r ( N + !) + T V(N + l) 2 + 2r(N + l) + l^j 
ft := ft(r) = y ^~ t 2 ( X + T (^ + !) + V(N + 1) 2 + 2t(N + 1) + 1 ) > 4 
so that, P_ 2 fir) < (3 2 fir )• 

We easily check that: (ft, ft) satisfies (12.481) and (12.491) if and only if, 

ft € (P v Pi) and ft = 2 + rft + y/( 2 + rft) 2 -ft (ft - 4(JV + 1)) := 99 + (ft), 
or equivalently: 


ft € (ft^ft) and Pi = 2(iV + l)+Tft + y/(2(N + 1) + rft ) 2 - ft (ft - 4) := ^(ft). (2.57) 
We prove the following: 

Theorem 2.12. // ^ / r £ ('7o>' r i) and we let pffir), i = 1,2, defined in (|4.54J) and (14.551) 
below, then the following hold: 

max{4(iV + l),^ 1 (r)} < Pfifir) < Pi fir) < ft(r), 
max{4 ,P 2 (t)} < Pfi fir) < Pi fir) < ft(r); 

moreover problem ( P) T admits a radially symmetric solution if and only if the pair (ft, ft) 
satisfies: 

Pi € (ft" (r),Pfi (t)) and p 2 = (pf (£ 1 ), (2.59) 

or equivalently: 

ft € {P 2 fir), Pt fir)) and ft = 9 ft (ft). (2.60) 

Furthermore , 

Pffir) ->• 4 (N + 2) as r -9 (2-61) 

□ 

We refer to Theorem 14. 121 Theorem 14.161 and Theorem 14. 22 1 for more detailed statements and 
a better grasp about the origin of the value: ft* 1 , i = 1,2. 

To establish the results stated above, we use a blow up analysis. 
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More precisely, for every «£lwe let (v\(r, a), V2(r, a)) the solution of the Cauchy problem 
for the radial ODE system corresponding to (P) T with initial conditions: 

ui(0, a) = a, ^(0, a) = 0, ui(0, a) = 1 ) 2 ( 0 , a) = 0. (2.62) 

We check that such solution is unique and globally defined for every r > 0, and it satisfies 
the required integrability conditions: \x\ 2N e vi G L 1 (M 2 ) and e V2 G L 1 (M 2 ). 

Thus by setting: 

/»+OO /*+oo 

f3i(a) := / r 2N+1 e V1 ^ dr and @2(01):= / re V2 ^ dr (2.63) 

Jo Jo 

we can use the scale invariance (|2.27l) . in order to see that the whole set of pairs (fJi,^) for 
which (-P)r admits a radial solution is fully described by the smooth curve: 

(^i(a),/3 2 (a)),a G M. (2.64) 

By using blow-up techniques introduced in [6]]!y58] in the context of Liouville-type equations, 
we shall be able to identify the limit value of /?*(«), as a —>• +00 and a —>• —00 for z = 1,2, 
and we check that indeed they yield to the (sufficient) statement about the existence of radial 
solutions claimed in Theorem 12.11 Theorem 12.31 and Theorem 12.91 The ’’necessary” part is 
contained in Theorem 14.221 

It is worth to notice that the ’’limit” values of /3i(a ) in general do not capture the sharp bounds 
for the pairs {(i\ . fa)-, and this fact has already come up for a single equation Liouville-type , 
as discussed in [59], and [25]. So our classification result gives also indication of certain non¬ 
degeneracy properties enjoyed by radial solutions of problem (P) r , to be further investigated. 

3 Some Useful Facts 

In this section, we collect some general informations about solutions of (I2.20p . (I2.2ip . in 
terms of the coupling matrix K and the ”flux-pair” (/3i, /? 2) 5 of independent interest. They 
generalise some known facts concerning the symmetric non degenerate case, see P21EE3H2I]- 
As above, we let: 


Ui(x) = Vi(x) + 2Ajlog |x|, * = 1,2 (3.1) 

and formulate (12.201) . (12.211) in terms of (vi,V 2 ) (the regular parts of ( 111 , 112 )), as follows: 

—Aui = k\ t i\x\ 2 Nl e V1 + ki^\x\ 2 N 2 e V2 in M 2 

< —Av2 = ^2,1 \x\ 2 Nl e Vl + k2 t 2\x\ 2 N 2 e V2 in M 2 

R — — f R„ — _L I |/v.|2Af2 p i>2 

— 27 T Jr 2 e > P2 — 2n Jr 2 l X l 6 
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(3.2) 















The following holds: 

Lemma 3.1. If (v i,v 2 ) satisfies (13.211 then we have: 


vfix) = log |x| + 0(1), for \x\ > 1 

rd r Vi(x ) —> dgvfix) —>• 0 as r = |x| —>• +00 


(3.3) 


with /3f° in (12.301) . i = 1,2; and (r,6) polar coordinates. In particular the condition: 


fii° — + kijfij > 2(iVj + 1), \/i j £ {1,2} 


(3.4) 


is necessary for the solvability of m- 
Furthermore, 

i) if kj j < 0 for i 7^ j £ {1, 2} then necessarily the coupling matrix K must satisfy (12.381) . 

ii) If kj.j > 0 for i / j £ {1, 2} but h h i < 0 for some i £ {1, 2}, then necessarily detK < 0. 

Proof. : By an obvious modification of the arguments provided in PS El, one can derive 
that every solution (vi,v 2 ) of (|3.2D must satisfy: vf £ L°°(R 2 ), i = 1,2. Hence, can use 
such information and suitable potential estimates, as in EQ3QZ] or in m, in order to deduce 
(13.31) . At this point, we can obtain (13.41) as an easy consequence of (13.31) and the integrability 
of the function \x\ 2Ni e Vi in M 2 , i = 1 , 2 . 

Concerning i), we observe that if kij = 0 for some i j £ {1,2} then (12.381) follows directly 
from (|3.4|) . While if k l: j < 0, Vi 7 ^ j £ { 1 , 2 }, then from (13.41) we see that k t _ t > 0 for all 
i £ {1,2}. Furthermore, we can rewrite (13.41) equivalently as follows: 


(&1,1&2,2 — &l,2&2,l)/3l > 2(A r i + 1)^2,2 + 2(iV2 + l)|fcl,2 
(&1,1^2,2 — fcl,2^2, 1)^2 > 2(A^2 + l)fel, 1 + 2(Aq + 1)1^2,l| 


(3.5) 


and (12.381) readily follows. Similarly, to deduce ii) we assume for example that k\ t 1 < 0, then 
from (13.41) we find that, 


fi 2 (detK) = @ 2 (k 2 , 2 1^ 1,1 1 + k 1}2 k 2 ,i) > 2(Ni + l)fc 2 ,i + 2(N 2 + l)|fci,i| > 0, (3.6) 


and necessarily detK < 0 as claimed. 


□ 


Lemma 3.2. Let K = (&ij)jj=i ,2 satisfy: 


k\2 ■ k 2 1 > 0. 

If problem ()2.21) admits a solution then the pairs (/3 i,/ 3 2 ) must satisfy: 


ki,i\k2,i\fii + k 2)2 1^1,2 \02 + 2 ^ 1,2 |^2,i|/3 i/52 — 4(lVi + 1 )|/c2,i|/3i — 4(A^2 + 1)1^1, 2 1^2 — 0. (3.7) 
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Proof. : In case k\ 2 = 0 = k 2 _\ then the left hand side of (13.71) is identically zero and we have 
nothing to prove. On the other hand, the system in this case decouples into two singular 
Liouville equations, whose solvability requires ka > 0 for every i = 1,2 and whose solution 
set is completely described in m- 

Furthermore if k \_ 2 = 0 but k 2 .\ / 0 (or fc 2) i = 0 but k\ t 2 / 0) then only the hrst (or second) 
equation becomes a singular Liouville equation and (13.71) simply leads to the following well 
known fact, 


ki,i/3i — 4(IVi + 1) = 0 (or k 2 ^ 2 f5 2 — 4(1V 2 + 1) = 0) 

see [IHISI]. 

So the true interesting situation to analyse occurs when, 


ki,2 • ^ 2,1 > 0. (3.8) 

Claim 1: If k\^ = k 2) \ / 0 then (13.71) holds, or equivalently: 

fci.i/?! 2 + k 2 rfl + - 4 (IVi + l)/3i - 4(iV 2 + l)/? 2 = 0. (3.9) 

It is clear that (E2D and (13.91) are equivalent when ki >2 = & 2 ,i 7 ^ 0. 

To establish (13.91) we use Pohozaev’ s trick and multiply the first equations in (13.21) by 
Vvi(x) • x, and the second equation in (13.21) by Vu 2 (x) ■ x, and integrate over the ball B r = 
{x € M 2 : \x\ < r} to obtain the following identity: 


L r (5 iv ” >|2 - (§£)) ^ = ^ L r2N ' +leV ' da 

- 2 (Ni + 1 )k ii i [ \x\ 2Ni e Vi + k i}j [ \x\ 2Nj e Vj Vvi(x) ■ x, i ^ j £ {1,2} 
J B r J B r 

where we have used the well known identity: 


(3.10) 


(Av(x))'Vv(x) 


x = div(yv{x){x ■ v{x) 



) 


and Green-Gauss theorem. 

On the other hand, if we multiply the hrst equation in (13.21) by Vu 2 {x) ■ x, and the second 
equation in (13.21) by Vui(x) • x, and then we integrate over B r , we find: 


— / A Vi(Vvj ■ x)dx = ki t i / \x\ 2Nl e Vi (Vvj(x) ■ x)dx 

J dB r J B r 

+ kij f r 2Nj+1 e Vj da — 2(Nj + l)kij f \x\ 2Nj e Vj , 1/j €{1,2} 

J 8B r J B r 


(3.11) 
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where again we have used the Green-Gauss theorem. 

Since k\^ = &2,i, from (13.101) and (13.111) we derive the following identities: 


k\ t 2 / (Avi(x)'VV 2 (x) ■ X + Av 2 (x)'Vvi(x) ■ x)dx = 


— &i, 2 &i,i [ \x\ 2Nl e vl ( x \VV 2 (x) ■ x)dx — k\^k 2,2 [ \x\ 2N2 e 02 ^ {S7 v\(x) ■ x)dx 
J B r J B r 

-k\ 2 f r 2Nl+1 e vl da — kf 2 [ r 2N2+1 e V2 da 

’ JdB r ’ JdB r 

+ kf 2 (%(Ni + 1) / \x\ 2Nl e Vl dx + 2(N 2 + l) [ \x\ 2N2 e V2 ^ x) dx] 

\ J B r J B r J 


and 


(3.12) 


* w L/fi |Vvi|a- (^) ) dc,+ki ’ i L r (^ v2]2 -(^) l*- 


- ki t ik 2 ,2 


IdBr 


r 2Ni+l e vi d(J + f r 2N 2 +l e v 2d(J 
JdBr 


—2(N\ + 1) [ \x\ 2 Nl e vl dx — 2(N 2 + 1) [ \x\ 2 N 2 e V2 ^dx 
J B r J B r 

+ k\ y 2k \,i f \x\ 2 Nl e vl ( x \'Vv2(x) ■ x)dx + ki } 2k2,2 [ \x\ 2 N ‘ 2 e V 2 ^ x \'Vvi(x) ■ x)dx 
J B r J B r 

By recalling the identity: 


(3.13) 


Avi(x)(Vv\(x) ■ x) + Av2(x)(Vv2(x) ■ x) = div(x ■ VuiVu2 — x ■ V ± vi'V ± V2) = 
= div(x ■ X 7 v 2VU1 — x ■ V^^V^ui), 


(3.14) 


where: we can turn also the right hand side of (|3.12l) into a boundary 

integral, and by summing up (13.121) and (13.131) we arrive at the following: 


ku L r (5 1 Vt ’ 2 ' 2 - (£)) d<7 +122 L r (5 iv " 112 - (£r 

+ki > L r {%% ~ (v± ” 1 ■ ,/)<v±! ’ 2 ■ "0 = 


\ x \ 2N 2 e V2(x) 


da 


(detK) f f r 2Nl+1 e Vl da + [ r 2N2+l e v2 da 
\JdBr JdBr 

- (detK) ( 2 (JVi + 1) [ \x\ 2Nl e vl + 2 (N 2 + 1) [ 
\ J B r Jh 


(3.15) 
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At this point, by taking into account the asymptotic behaviour of v\ and v 2 as given in (13.311 
and (13.41) . we can pass to the limit, as r —>• oo in (13.151) . and conclude: 


K 2 PTP 2 - \h 2 ^T) 2 ~ = —{detK)(2(Ni + l)/?i) + 2 (N 2 + l)p 2 ) (3.16) 

with, 


PT = h,iPi + k h2 p 2 and /3f° = k h2 Pi + k 2>2 p 2 , (3.17) 

see (12.301) . Identity (13.161) reaffirms i) and ii) of Lemma 13.11 

By inserting (13.171) into (13.161) and by carrying out straightforward calculations, we easily 
arrive at provided that: detK / 0. 

On the other hand, by Lemma f3.ll we also know that detK 7 ^ 0 enters as a necessary condition 
for the solvability of (13.21) . except in the cooperative case (I2.26P where K could be degenerate. 
But for degenerate cooperative systems, the matrix K satisfies: 


ki,j > 0 i, j E { 1 , 2 } and fei,ife 2 ,2 = kj 2l 


and so necessarily, 


000 k\ 2 QQQ 1 k\ 2 

P 2 = t^Pi and v 2 = - — v x . 

fepi fei,i 


In other words, by setting: 


P = P?° v = vi and a = 

«i,i 

we see that in this case (13.21) reduces to the following single Liouville type equation: 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


—Av = £q, i\x\ 2Nl e v + ki :2 \x\ 2N2 e av in R 2 
P = ^ fu2(k i,\\x\ 2Nl e v + k lj2 \x\ 2N2 e av )dx 
with feu > 0 , fei j2 > 0 and fegife 2i2 = fe 2 2 - 

Equations of this type arise in the construction of Self-gravitating Cosmic Strings (cfr. [TT]h 
and have been analysed in para Eli eh]. 

In particular, in this case, for i = 1 we may complete (13.101) . by using (13.191) and (13.201) and 
obtain that every solutions of (13.211) satisfies: 


) da=k '- i L 

— 2(N\ + l)fei ; i / 

Je 


r^+l e v d(T + h^ 


r 2N 2 +l e av da 


a JdB r 


x \ 2Nl e v -k l2 2{N2 + l) [ 
a Je 


\x\ 2N2 e av . 


B r 


22 

























Therefore by passing to the limit in the above identity as r —» oo, we conclude the following 
identity: 


P 2 = 4(JVi + l)k lt 1 P 1 + 


4 (N 2 + 1) 


kl.2^2 


(3.22) 


with 


1 

/ 

\x\ 2Nl e v -- 

1 

f \x\ 2Nl e Vl , 

“ 2tt , 

Jr 2 


2 vr j 

R 2 

1 

f 

\x\ 2N2 e av 

1 

[ |z| 2A V 2 

“ 2tt , 

Jr 2 


“ 27T 

Jr 2 


(3.23) 


P = h,iPi + h,2p2- 


At this point, by recalling that a = and fci,ifc 2j2 = k\ 2 i we easily derive (13.91) for P\ and 
P 2 simply by inserting (13.231) into (13.221) . 

Thus identity (13.91) is established in all cases, and the proof of the claim is completed. 

To conclude the proof we shall see how to use (13.81) in order to reformulate the (non- 
symmetric) system (13.21) into a symmetric one, for which we can use (13.91) . To this purpose, 
we set, 


Di(.x) = v\(x) + log 1 ^ 2 , 11 V 2 (x) = v 2 (x) + log |&i )2 | 

and observe that, if (vi,v 2 ) satisfies (13.21) then (v\ , v 2 ) satisfies a similar problem with the 
symmetric coupling matrix: 


( £11 +i 

K=( l fc -l f 

\ +1 fc 22 

V IM 


(3.24) 


and where the ± sign is chosen according to the sign of fci i2 (or equivalently k 2j \). Further¬ 
more, 


Pi : = 


1 


x\ 2Nl e Vl = \k 2 ,i\Pi and p 2 := 


-! 

J r: 


x\ 2N2 e V2 = |fei j2 |/3 2 . 


27T J^2 1 27T 

Hence, we can apply the Claim to (Pi, $ 2 )- Thus, by virtue of (13.241) . from (13.91) we find: 


(3.25) 


^rPl + ± P 1 P 2 ~ mi + l)3i - 4(iV 2 + 1 )P 2 = 0, 

F2.ll |«1, 2 | 

and so from (13.251) and (13.261) . we readily derive (13.71) . 


(3.26) 


□ 
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Remark 3.3. It is easy to check that, in the symmetric case, i.e. k\- 2 = k 2 ,i, the following 
identity holds: 


(ki,ik 2 , 2 -k 2 12 )[(N 1 + l)p 1 + {N 2 + l)p 2 )} = 

’ (3.27) 

[(iVi + 1 )k 2 , 2 - ( N 2 + l)fci, 2 ]/?r + i(N 2 + 1 )&!,! - (N! + 1 )k lt2 ]P?, 

and we can use it together with (13.161) to deduce the following equivalent formulation of (13.91) 
expressed in terms of (3f°, i = 1,2; 


k 2 , 2 (/3n 2 + fci,i(^°) 2 - 2fci, 2 /3 1 00 /3 2 00 - 4[(A r 1 + l)fc 2>2 - (N 2 + l)k lt2 ]fi° 

(3.28) 

- 4[(1V 2 + l)fci,i - (Ni + l)fci, 2 ]^° = 0, 

which has the advantage to holds even for a degenerate coupling matrix K. 

Arguing as above, one can derive a similar identity for the non-symmetric case, provided that 
k\- 2 ■ k 2: i > 0, we omit the details. 

Clearly, Lemma 13.1 1 and f3. 21 imply Proposition 12.41 and Corollary 12.61 as stated in the previous 
section. 


Next, we turn to analyse the radial solvability of (13.21) . To this purpose, we observe that a 
radial solution (about the origin) (ui(r), v 2 (r)) of (13.21) . may be characterised equivalently as 
satisfying the following boundary value problem: 


— (rv'^r))' = Aq,i r 2Nl+1 e vl ^ r ' > + ki t2 r 2N2+l e V2 ^ r \ r > 0 
—(rv' 2 (r)y = k 2 ir 2N i+ 1 e v i ( r ) _)_ 2 r 2N2+1 e V2 ^ r \ r > 0 

< 

u'(0) = 0, lirn^oo ru'(r) = 

Vi € C([0,oo)), i = 1,2, 
with /3f° i = 1,2, defined in (12.301) . 


To this purpose, we consider the Initial Value Problem associated to the system of ODE’s in 
(13.29[) . and show that it admits a unique globally defined solutions which also accommodate 
the boundary condition required in (13.291) . To be more precise we let, 

fi(r)= [" s 2Ni+1 e Vi ^ s) ds (3.30) 

Jo 

and observe that we can express any local solution of the system of ODE’s in ()3.29|) defined 
in the interval I = [0,72), equivalently as follows: 
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rv[(r) + fci,i/i(r) + 2 / 2 ( 7 -) = 0, 0 < r < R, 

< rv' 2 {r) + k 2 , 1/1 (r) + k 2j2 f 2 (r) = 0 < r < R, (3.31) 

y i (0)=0,v i eC([0,R)), i = l,2. 

We have: 

Lemma 3.4. If the coupling matrix K is symmetric and (vi(r),v 2 (r)) satisfies (13.311) with 
fi = fi ( r ) in (13.301) . then: 

l) r 2(W+ l) e vi(r) + r 2{N 2 +l) e v 2 {r) _ 2^ + l)f^ r ) - 2(N 2 + 1 )/ 2 (j") + 

i (3.32) 

^(ki,ifi(r) + 2fei,2/i(r)/ 2 (r) + A: 2)2 /|(r)) = 0, V0 < r < 

ii) _ 2 (W + l)/^(r) + (^//(r) = - [ kijfi(s)fj(s) ds, 

2 ./o (3.33) 

VO < r < R, Vi j e { 1, 2}. 

Proof. The above identities were pointed out in |59| for more general systems, and are ob¬ 
tained simply by computing the derivative of the term in the left hand side of (13.321) and 
(13.331) respectively, and by using (13.311) . 

□ 

Again we observe that, as in the proof of Lemma 13.21 it is possible to derive identities anal¬ 
ogous to (I3.32p and (13.331) . also in case we drop the symmetric assumption on the matrix K, 
but we assume instead that: k\^ ■ k 2) i > 0, we omit the details. 

Furthermore, for solutions of (13.291) . we can deduce (12.241) and (|2.25l) as a direct consequence 
of 113.321) and (13.331) . just by letting r —»• oo. 

Remark 3.5. In case K is strictly positive definite, (a necessary condition when k\, 2 = 
k 2) i < 0), from (13.321) we obtain that the monotone functions /i(r) and f 2 (r ) in (|3.30l) are 
uniformly bounded by a constant depending only on N\ and N 2 . Therefore, by a standard 
Picard iterative scheme, we see that any (local) solution of the Cauchy Problem associated 
to the system of ODE’s in (13.291) can be globally extended in [0,oo) to define a solution of 
(13.291) with fii := lim^oo ffir), and (3f° defined by (12.301) . i = 1,2. 

4 Main Results and their Proof 

We devote this section to analyse radial solutions for problem ( P) T with r e (0,1), or 
equivalently to investigate solutions (ui(r), v 2 (r)) of the problem: 
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— (rv[(r))' = r 2N + 1 e v i ( r ) _ rre V2 ^ r \ r > 0 , 

— (rv'r,(r)Y = re° 2 ^ — Tr 2N+i vi(r) > q 

< 2V " ’ (4.1) 

v'(0) = 0, lim^oo ru'(r) = -0j, r , 

Vi € C([0,cx))), i = 1,2, 

with 

Pi,T’= Pi-Tfo and 0 2 , r ■= 02 ~ t0i, (4.2) 


and 


or equivalently, 


/*00 7*00 

01 := / r 2 iV+ 1 e t ’ 1 (ir and 0 2 := / re V2 dr\ 

Jo Jo 


ft = %±I&1 and & = fc + rft,, 

1 "T" 


(4.3) 


(4.4) 


1 — r 2 ' 1 — r 2 

According to the results established in the previous section, we know that the following 
conditions are necessary for the solvability of (14.11) : 


0i 2 + $ - 2t/3i/3 2 - 4(AT + l)ft - 4/3 2 = 0, (4.5) 

01 > 4(1V + 1), 0 2 >4, (4.6) 

01,r =01-T0 2 >2(JV + 1), 02,r = 02 - T01 > 2 ; (4.7) 


and we recall that (US) can be expressed in terms of 0i iT and 02 , T equivalently as follows: 


01,r + 02, r — 2'T01,r02,r — 4[(iV + 1 + t)0 i, t + (1 + t( 1V + l))02,r], 


(4.8) 


or as follows: 

01 (01,r - 2(iV + 1)) + 02(02,r - 2) = 2(A^ + 1)0 X + 20 2 . (4.9) 

Our main goal in this section is to investigate to what extent the conditions (14.5[) . (14.61) and 
(14.71) are also sufficient for the solvability of (14.11) . 

To this purpose we recall that problem (14.11) is invariant under the following scaling property: 


v\ (r) —»■ vi t \(r) = vi(Xr) + 2(N + 1) log A 
v 2 (r) —)• v 2y \ (r) = v 2 (\r) + 2 log A, 


(4.10) 


in the sense that: (vi,v 2 ) solves ( P) T if and only if (u^a, u 2 .a) solves ( P ) T . We shall distin¬ 
guish between this 1 -parameter family of solutions by their initial value at r = 0 . 
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Notice also that, for r E (0,1), the identity (13.321) for problem (14.11) reads as follows: 


(Jo + (fj teMt) ~ 2r (fj t 2N+1 e V1 ^ dtj ^jT te V2 ^ dt 

-4(N + 1) t 2N+1 e v lW dt) - 4 f jT te V2 ® dt) + 2r 2 ( 7V+ V l(r) + 2 rV 2(r) = 0, 


(4.11) 


for all r > 0 . 


From (14.111) we obtain the following uniform estimate about solution of (14.11) . independently 
from their initial data and the values of (3\ and 02 ■ 

Lemma 4.1. Let r E (0,1) and (rq, r^) be a solution of (14.11) . There exists a suitable constant 
C = C(t,N) > 0, depending only on r and N, such that: 


v\ (r) + 2 (N + l)logr < C and V 2 (r) + 21ogr < C, for any r > 0. 

Proof. By direct calculation, from (14.111) we easily derive the following estimate: 

r 2 {N+l) e vPr) +r 2 e v 2 {r) < 2[(1V + l ) 2 + 2t(N + 1) + 1] ^ 

1 — T 2 

and by (14.131) . we readily derive (|4.12[) . 

An important and useful consequence of (14.121) is given by the following: 


(4.12) 


(4.13) 

□ 


Lemma 4.2 (Harnack’s inequality). Let r E (0,1) and (v\,V 2 ) be a solution of (14.11) . For 

Rq 

every 0 < ro < Ro, there exists a constant Co = Cq(t,N, —) > 0 (depending only on t, N 

ro 

R() / R ^ \ 

and — ) and 7 = 7 ( — ) E (0,1) such that, 
ro \ r o J 


max v\ < 7 min v\ + 2 (N + 1)(7 — 1) log ro + Co, 
[ro,J?o] [u>,-Ro] 

max V 2 < 7 min V 2 + 2(7 — 1 ) log ro + Co- 
[r 0 ,i?o] [ro,Ro] 


(4.14) 


Proof. Set Lq : = 


Rq 

ro 


and define D 


2 <r< 



We rescale ( 77 , tq) as follows: 


v\(r) = ri(ror) + 2(N + l)logro and V 2 (r) = V 2 (r$r) + 21ogro, (4-15) 

so that, in view of (14.101) . the pair (hi, ^ 2 ) still satisfies (14.11) together with (14.121) . In particular 
by setting: 


/i(r) = r 2N+1 e Vl P ) — rre V2 and hi r ) = re v — rr 2N+1 e Vl , 
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by (14.121) we find a constant C = C(t, N ) > 0, depending only on r and N, such that: Vi < C 


in D, and 


fi 


<C,i = 1 , 2 . 

L°°(D) 

Therefore, for i = 1,2, we let ifi the unique solution of the following Dirichlet problem: 


-(rip'J = fi in D 

ipi = 0 on dD 


* = 1 , 2 , 


(4.16) 


and derive, by standard elliptic estimates (see [3D]), that, 

fi 


max |^j| < A 
D 


L°°(D) 


with a suitable constant A = A(Lq) >0. As a consequence, or i = 1 , 2 , the function 
4>i := ifji — (vi — C ) defines a positive harmonic function in D , for which we can use Harnack 
inequality in D' :={1 < r < Lo} CC D, and obtain: 

sup 4 >i < — inf (j)i, i = 1 , 2 . 

D' 7 D ' 

with a (universal) constant 7 = 7 (^ 0 ) € (0,1). In other words, 

sup Vi < 7 inf Vi + (1 — 7 )C + (1 + 7 ) max |V>i| 

O' D' 

with 7 = 7 (^ 0 ) £ (0,1), and the desired estimate (14.141) follows. 

□ 

Finally, we observe that, in view or Remark 13.51 for r € (0,1) and for any ( 07 , 02 ) E M 2 , the 
initial value problem: 


— (rv'])' = 7 ’ 2 iV + 1 e 1, i — rre V2 

— (rv oV = re 1 ' 2 — rr 2Ar + 1 e ,; i 


(4.17) 


Ui(0) = Oj, u'(0) = 0 * = 1,2 

admits a unique solution defined for all r > 0 , and such that: 


r 

Jo 


r c - 

r 2 N+ 1 e vl dr < C and / re V2 dr < C 

Jo 


with a suitable constant C = C(N,t) > 0 depending only on N and r, (but independent of 
the initial data ( 01 , 02 )). 

Remark 4.3. By virtue of the scale invariant property (14.101) . it follows that the set of pairs 
(/3i, /5 2 ) for which (14.11) . (14.21) admits a solution is fully described by the curve: 


(/3i(o),/3 2 (a)), ol E 


(4.18) 
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where, 


POO POO 

Pi{a) = / r 2N+1 e vl ^ a) ds, p 2 {a) = / re V2 ^ a) ds, 

Jo Jo 


(4.19) 


and a), v\ (r, a)) is the unique solution of the initial value problem (|4.17l) with initial 

data specified as follows: 

a\ = a and a 2 = 0. (4.20) 

Our main purpose will be to identify the limiting values of /3i(a), i = 1, 2, as a —>• ±oo. 

To this end, we start to give a detailed description about the portion £ of the ellipse defined 
by (13.321) which is contained in the first quadrant of the (pi,p 2 )-Y>\ ane, namely: 


£ 


fil + fil - 2 rp 1 p 2 - 4 (N + 1 )p l - 4/5 2 = 0 
p l >0,p 2 >0. 


(4.21) 


For this purpose, it is useful to identify first the intersection of £ with the line: Pi — tP 2 = 
2(N + 1) and the line: Pi — tPi = 2. To this end, we set: 


(.N + l) 2 + 2r(lV + 1) + 1 =: D(t, N ) 
and by direct calculations, we verify the following: 


(4.22) 


Pi>o,p 2 >o 

< Pl + Pl - 2rpip 2 - 4 (N + 1 )Pi - 4/3 2 = 0 «=> < 

Pi - rp 2 = 2(N + 1) 

and similarly: 


Y~~f2 ( N + 1 + r + = : §_i( T ) 

j r 2 (l + r (^ + 1) + y/~D^J =: P 2 (t) ; 

(4.23) 


Pi >o,P 2 >0 

< pl + Pl - 2rpip 2 - 4(N + l)Pi -4p 2 = 0 
P 2 -tPi = 2 

Furthermore, we can explicitly solve the quadratic equation 
this way, we are lead to consider the functions: 



1 _ t2 [n + 1 + r + yfE) —: Pi(t) 
j^ r 2 ( 1 + r ( iV + !) + t'/d) =: P 2 (t). 

(4.24) 

in (13.321) in terms of Pi or p 2 . In 


ViiPi) = 2 + tPi± \J{2 + tPi) 2 - Pi (Pi - 4(N + 1)), with 0 < Pi < Pfir) (4.25) 


and 
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<pf(P 2 ) = 2(iY + l)+T/3 2 ±\/(2(iV + 1) + t/3 2 ) 2 - P 2 (P 2 - 4), with 0 < /3 2 < /3 2 (t). (4.26) 


To simplify notations, from now on we shall drop the dependence of P-(t) and Pfi r) i = 1 , 2, 
on the parameter r. 

By direct inspection, one can easily check the following: 

Lemma 4.4. We have: 


(*) <Pi 




is strictly monotone decreasing with inverse 


(vt) 1 — vt '■ \§_2^P'A \P_i^Pi]'i 


(**) Ti : [4(AC + 1), /5i] ~^ 0) P_ 2 is strictly monotone increasing; 


(Hi) ip 2 : [4 ,p 2 ] ->• 0,^ 


is strictly monotone increasing. 


□ 


For later use, we emphasise that, 


<Pi(P l) = P 2 ’ 

¥>2 (A 2 ) = £i> 

and moreover, for (/3i, /3 2 ) £ £ there holds: 


(4.27) 

(4.28) 


0 < Pi < Pi and 0<p 2 <p2- (4.29) 

Next, by recalling (14.91) . we observe that (Pi tT — 2(N + 1)) and (/3 2 , T — 2) cannot be simulta¬ 
neously negative. Thus we obtain the following description of £: 

Lemma 4.5. The pair (Pi, P 2 ) € £ if and only if it satisfies one of the following set of 
conditions: 

(' i ) Pi,r > 2(iV + 1) and p 2 , r >2, Pi € (p_ v Pi) and p 2 = Ti(Pi) 

V ' , (4-30) 

(or equivalently, p 2 E \P_ 2 -,P 2 ) an d Pi = T 2 ($ 2 )); 

(H) Pi, T < 2(AC + 1) and p 2 , r > 2, p 2 E (4,/3 2 ] and Pi = tpf (p 2 ) E (0,ft]; (4.31) 


(in) Pi >T > 2(N + 1) and fi 2 , r < 2, ft £ (4(A7 +1), /?-,_] and p 2 = ip x (Pi) E (0,/^]. (4.32) 
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Proof. From (14.911 we see that (14.30I) - (I4.32I) cover all the possibilities, and the relative expres¬ 
sion for (Pi , fi 2 ) can be easily verified by virtue of the definition of and ip 2 . 

□ 

Observe that, ( 14 . 301 ) accounts for the conditions ( 14.511 and ( 14 . 71 ) . In order to account also for 
m, we search for the intersections of the set £ in ( 14.2111 with the lines p\ = 4(IV + 1 ) and 
/?2 = 4 respectively, and find: 


Pi + pi - 2rPip2 - 4 (N + 1 )p 1 - 4/3 2 = 0 
Pi = 4 (IV + 1 ) and /3 2 > 0 


P‘2 


4 + 8 t(N + 1) 


PUr) (4.33) 


and 


Pi + Pi - 2 rpip 2 - 4 (N + 1 )Pi - 4/3 2 = 0 
P 2 = 4 and p 2 > 0 
Also, we dehne: 


Pi =4(IV + 1) + 8 r =-.p{ (r). (4.34) 


PT (r) ■■= 2 rp* 2 (r) = 8 r (1 + 2t(N + 1)), P* 2 * (r) := 2 T p{ (r) = 8 r (N + 1 + 2r) , (4.35) 

and observe that, not only (4(A r + 1 ),/ 32 (t)) and (Pl(r), 4) € £, but also: 

(Pi*(r),P* 2 (r)) and (^(r), P* 2 * (r)) € £. (4.36) 

We expect such values to play a role in identifying the pairs (/3 i,/ 3 2 ) for which (14. 1 II - (14.311 
admits a solution. 

Proposition 4.6. We have: 

(i) 


PT(t) = 4 (N + 1) « ^(r) = 4(JV + 1) and p* 2 (r) = P 2 (t) 


N + 1 


T = 


/n particular, 


1 + y/1 + 4(N + 1) 


/3r (r) - tP* 2 (t) = 2(N + 1)^t = . 


- r (1) 
2 — T o • 


(**) 


/3**(r) =4/3 (r) = 4, and /3*(r) = /^(r) 


-2 

r = 


1 


_. A 2 ) 


N + 1 + ^/ (N + l ) 2 + 4 


= : t. 


In particular, 


P 2 * ( T ) ~ T P* ( T ) = 2 <=>■ r = Tq 2) . 


(4.37) 


(4.38) 


(4.39) 


(4.40) 
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Moreover: 


(4.41) 


-( 2 ) 


0 < 7q < Tn < 


( 1 ) 


0 


Proof. By straightforward calculations we see that, ft*(r) = 4(1V +1) if and only if r = 
Furthermore recalling that: ft*(r) = 2rftj (r), we also find that 4(1V +1) — r/31(r) = 2 (IV +1) 
if and only if r = Tq X \ and therefore: ftft) = 4 (IV + 1) and ft (r) = /3 2 ('r) if and only if 


t = T, 


(i) 


So (z) is established, and (zz) follows exactly in the same way. 
Finally to check (14.411) . we observe that the function f(t) = 


for t > 0. Hence, for N > 0, we have: 


t + Vt 2 + 4 


, is strictly decreasing 


.( 2 ) 


= /(JV + l)</ 


1 


JV + 1 


= T, 


(i) 


□ 


As an immediate consequence of Proposition 14.61 we obtain: 

Corollary 4.7. For every t G (0,1) we /zane: 

— ( 2 ) 

(z) 0 < ft(r) < ft(r) and equality holds if and only if t = Tq . 

(zz) 0 < ft(r) < / 3 2 (t) and equality holds if and only if r = Tq 1 ^. 

Proof. In view of (14.291) and (14.271) . we see that properties (z) and (zz) follow easily from 
Proposition 14.61 . 

□ 


(4.42) 


Corollary 4.8. There holds: 

(0 

ft (ft < 4(IV + 1) 4(IV + 1) - rft(r) > 2(IV + 1) > ft* (r) - rft(r) 

^ ft* (ft < 4(IV + 1) ^ r G (o, rj 1} ) . 

Moreover, for every r G ^0, Tq 1 ^ we have: 

ft (ft < 4(IV + 1) = ftj~(ft(ft) ( or equivalently ft(r) = ftft4(IV + 1)) < ft (ft) 
and ft*(r) = ft"(ft (ft) < ft (ft- 

(4.43) 


(zz) 


ft (ft < 4 «=► 4 - rft(r) > 2 > ft*(r) - rft(r) « ft*(r) < 4 

.( 2 ) 


r G 0 , r ( 


(4.44) 


32 











Moreover, for every r G ^0, Tq 2 ^ we have: 


P 2 (t) <4 = <pf (P\ (r)) (or equivalently P{ (r) = (4) < Pi(t)) 

and PI *(t) = ipf(PUT)) < P 2 (r). 


(4.45) 


Proof. We prove ( i ) since (ii) follows in a similar way. The properties in (14.4211 follow by 
Proposition 14.61 and the monotonicity of the functions involved. To show (14.431) . we observe 
first that, for Pi = 4(iV + 1) and P 2 = /3|(r) we know that, pi — rp 2 > 2(N + 1) and we 
can also check that, P 2 — r/?i > 2, Vr G ^0, r^ 1 ^. So, we can use Lemma [4.51 together 
with the monotonicity of the function ip(( and ip 2 given in Lemma 14.41 to conclude that, 
/3 1 (r) = P 2 (P 2 ( t )) < P 2 (P 2 ) = 4(iV + 1), as claimed. The dual statement in (14.431) follows 
similarly. Furthermore, since for Pi = P\* and P 2 = P 2 we know that, Pi — TP 2 < 2 (N + 1), 
Vr G ^0, Tq 1 ^, and at the same time we can check that, P 2 — r/3i > 2, so we can use again 
Lemma 14.51 to derive: /3** = ip 2 (P 2 ) < /3 , as claimed. 


□ 


/ 2 "\ 

Remark 4.9. From the above discussion, we see that for t G (0, Tq the (necessary) condi¬ 
tions 63|). 631) imply the integrability condition 63b as it happens for the cooperative case 

( 2 ) 

r < 0. Thus, we expect that, when r G (0 ,Tq '], then the same uniqueness and non-degeneracy 
properties (as established in Theorem 12.,91 for r < should remain valid for radial solutions 
0f(Pr). 

Next, we wish to describe the relations between P**(t) and P 2 (t) when r G and 


P 2 * (t) and PI(t ) when r G ^r^ 1 1 j • Since for r G > 1 j we can always guarantee that, 


PT(t)-tP* 2 (t)>2(N + 1), 


(4.46) 


thus, we need to investigate when we can also realise the second integrability condition: 


P* 2 (t)-tP* 1 *(t)> 2. (4.47) 

Similar considerations can be applied about the relations between P((t) and P 2 *(t). 
Proposition 4.10. (i) There exists a unique value r-j 1 * G (^^ 5 ) such that, 

PU T ) ~ rPT(r) = 2, that is: P 2 (t) = P_ 2 (t) and P{*(t) = /3 x (r) r = rj 1) . 

More precisely, 

/^2 ( r ) ~ T Pi* ( r ) > 2 <(=> r G (0, t} 1} ) 

and 

0 <PT(r)<-pi(r) Vr G (0,1) \ jrp)} . 
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(ii) There exists a unique value t[ 2 ^ E ( 4 , -^) such that P\(t) — t(I 2 *(t) = 2 (N + 1), that 
is: PI(t) = /^(r) and P 2 *(t) = /3 2 (t) r = t} 2) . 

More precisely, 

Pl(r) - tP?(t) > 2(1V + 1) ^ t E (0, t [ 2) ) 


and 


Moreover, 


0 <P?(t)<P 2 (t) Vre( 0 ,l)\{r 1 (2) }. 


1 


( 2 ) / ^( 1 ) 


- < T{' < TM < 


1 

2 ' 1 ' 1 " 71 ' 

Proof. We establish (i) with the help of the function, 

Mr) = \ (P *2 (r) - rft*(r)) - 1 = 2(1 - 2t 2 )(1 + 2 r(N + 1)) - 1. 

We readily check that, for r = ^ +~i)^ - e (0, 5 ), the function if\ is increasing 

in (0, r] and decreasing in (r, + 00 ). Since ifi(\) = N + 1 > 0 while = —1, we find a 

unique : ifi(r^) = 0, and moreover ifi(r) > 0, Vr E (0, r^)- 

Finally from (14.291) we know that, 0 < P**(r) < /3 1 (r) and (i) is established. 

We observe that (ii) follows in a similar way by considering the function 

Mr) = \(P*i(r) ~ t/3 2 "(t)) - 2 (N + 1) = 2(1 - 2 t 2 )(N + 1 + 2 t) - (N + 1), 


we omit the details. 

hus i 
1 1 


Thus it remains to show that: r[ 2 ^ < r| 1 ^. To this purpose, we easily check that, Vr E 


, . there holds: ^ 2 ( r ) < V ; i( T )- 

2 v2/ 

So the desired conclusion follows by recalling that Tf is the unique positive zero for if i(r), 


while r{ 2 ^ is the unique positive zero for the function: if 2 (r), and E 


1 1 

2’7! 


,* = 1 , 2 . 


□ 


Corollary 4.11. (i) If t E (777'^, ^en 

01*( r ) = vt (02( r )) or equivalently P 2 = Ti (Pi*( T )) ■ ( 4 - 48 ) 

In particular, 

A(N + 1) < P^(r) < Pi*(r) < Pi(t) and P 2 (t) < P 2 (t) < P 2 (t). (4.49) 
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(ii) If t E (tq 2 \ rj 2 '* j then 

nr = ipiifii) or equivalently ft* = (4.50) 

In particular, 

4 < PJj) < 02*( r ) < A-M and P 1 ( T ) < 0l( T ) < 0l( T )' ( 4 - 51 ) 

Proof. If t E (r 0 ( 1 ) ,Ti (1) ) then we see that (14.461) . (14.471) hold simultaneously, and so we 

deduce (14.481) as a consequence of (14.301) . Furthermore, as r > Tq 1 ^ then /3|(r) < /3 2 (t - ), and 
so: 0**(r) = PtiP l( r )) > ^(/ 3 2 ( t )) = /^(t). Similarly, as r < r x (1) then ^*(r) < /^(r) 
and so /3|(r) = ^(/3|*(r)) > ^(/3 1 (r)) = /? 2 (t), and also (I4.49P is established. Part (ii) 
follows exactly in the same way, with the obvious modifications, we omit the details. 

□ 

Before we discuss our existence results for problem ( P) T we make the following simple obser¬ 
vation: 

/f(r)<^)^re( 0, |) * = 1,2 (4.52) 

and in particular, 

#*0-) = 0*(r) = 4(iV + 2) « r = i (4.53) 

In other words, when r = \ and problem (P) T — i reduces to the 2 X 2 -Toda-system, then 
for i = 1 , 2 , the values /3** and /?*, coincide with the only value allowed by solvability, see 
(I23TD . 

Indeed, we show below that actually, the values (I** and fl* i = 1,2, capture in a crucial way 
the (radial) solvability for (P) T for any r E (0,1). 

To this purpose, and in account of Proposition \TM Proposition 14. 1UI and (14.52|) . for r € (0,1), 
we define: 


01 (t) = < 


4(IV + 1), 0 < r < r, 




ip 


PI (t), |<r<r 1 


< r < 2 
( 2 ) 


^i(r), 


r{ 2) < r < 1 


and similarly, 


4, 0 < 

r < r< 2 > 

nr(r), 

T o 2) < T < 1 

02 ( T )5 

4 < r < r x (1) 

0 2 ( 0 > 

r{ 1} < r < 1 


Pi( T ) 


0* (r), 0 <t< 4 

< 01 **(r), 4 <r<r 1 (1) 

k 0l( T )> 'Tl (1) <'T<l 


02 " CO 


02 (0> 0<T< 4 

< 0r(r), 4 <T<r 1 (2) 

^(O. ^ 2 ) <T <1 


(4.54) 


(4.55) 
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We see that, (r), i = 1, 2 is a continuous functions of r, and for t ^ we have: 


max{^ 1 (r),4(AT + 1)} < P 1 (r) < Pf(r) < P x (t) (4.56) 

max{/? 2 (r), 4} < pf (r) < /3^(r) < /3 2 (t), (4.57) 

while, 

if r = ^ then /3f Q = Pf = 4(JV + 2 ), z = 1, 2 . (4.58) 

We prove the following: 

Theorem 4.12. If r € (0,1) and r / 4 ; problem (P) T admits a radial solution if and only 
if the pair (Pi,P 2 ) satisfies the following: 

Pi € (Pf(r),Pf(T)) and /3 2 = (fit (Pi) (4.59) 

or equivalently, 

P 2 e (/3 2 _ (t), ^ 2 "( r )) and Pi = V 2 (# 2 ), ( 4 -60) 

with pf- defined in (14.54[) and (14.551) . z = 1,2. 

The proof of Theorem 14.121 will require several preliminary steps, but first let us point out 
some of its interesting consequences. Indeed, in view of the definition of fif(r), j = 1,2 in 
(I4.54p and (14.551) . we deduce the following: 

Corollary 4.13. If r E (0 ,Tq ; ] then the conditions (14.51) . (14.61) imply the integrability con¬ 
dition 621). and they are necessary and sufficient for the radial solvability of (P) T - 

Proof. Simply observe that in this case: Pf(f) = 4(IV+1) and pf (t) = 4, so the pairs (Pi, p 2 ) 
satisfying: fix € (Pf(r), Pf (r)), p 2 = ipf (Pi) (or equivalently p 2 € (Pf (r), P% (r)), Pi = 
(p ~2 (Pi) simply describe the portion of the ellipse £ in (14. 211) contained in the quadrant: Pi > 
4(N + 1) and p 2 > 4. □ 

As already pointed out in Remark 14.91 when r € (0, Tq ) we expect that problem (P) T 
should keep also the same uniqueness and non-degeneracy properties established in [59] for 
cooperative systems, and which apply here for r < 0. 

Corollary 4.14. If r € [ti , 1) then the conditions (14.51) . (14.71) imply the condition (|4.6I) 
and they are necessary and sufficient for the solvability of (P) T - 
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Proof. In this case, Theorem 14.121 ensures the existence of a (radial) solution for ( P) T 
for any pairs (/3 i,/ 3 2 ) such that j3\ € (fl (t), /3 1 (t)) and /3 2 = pf(/3 1 ) (or equivalently 
/3 2 £ (/3 9 (t), /3 2 {t)) and f3\ = y?2~(/3 2 )), which fulfil exactly the portion of £ in (l4.21l) sub¬ 
ject to the constraint (14.71) . Furthermore, by Corollary 14.81 we see that the condition (14.61) is 
automatically satisfied in this case. 

□ 

Remark 4.15. Since for t € [rf 1) the condition (14.61) is ensured automatically by (14.51) 
and ED which also guarantee the existence of a radial solution (see below), thus we suspect 
that, whenever solvable, problem ( P T ) should admits only radial solutions. 

To establish the existence of radial solutions as claimed in Theorem 14.121 we analyse the 
Cauchy problem: (14. 171) - (14.201) . To be more precise, for a € R we let (ui(r, ct), u 2 (r, ct)) be 
the unique (global) solution of the Cauchy problem: 

—(rv'i)' = r 2N+1 e Vl — rre V2 , 

— {rv^y = re v2 — T r 2N+1 e Vl , 

< 

ui(0) = a u^(0) = 0 

u 2 (0) = 0 Un(0) = 0 

and set, 

poo poo 

f}\(a) = / r 2N+1 e vl dr, /3 2 (a) = / re V2 dr. (4.62) 

Jo Jo 

Clearly, the pair (/3i(a),/3 2 (a)) satisfy (14.5p . (14.61) and (14.71) . and in particular, /3j(a) is 
uniformly bounded with respect to a € M, for i = 1 , 2 . 

On the basis of Proposition 14.61 Proposition 14.101 and their consequences, we define: 


r > 0 
r > 0 


(4.61) 


and, 


a , x )P*( T ) for r G (0,r} 2) ) 

Pi,-oo(r) = < , . 

[^(r) for r £ [t{ , 1) 


# 2,-00 (r) = tpf (/3i,-oo (t)) 
similarly, we let: 


4 for r € (0, Tq 2) ] 

< P?(t) for Tetr'Vf) 
J3 2 (t) for r € [t{ 2) , 1) 


^(r) for r G (0, t| 1} ) 
^ 2 (r) for r € [t} 1) , 1) 


(4.63) 


(4.64) 


(4.65) 
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and 


0i,+oo( t) = V2 (02,+cxd(t)) 

It is important to observe that, 


4(IV + 1) for r € (0 ,Tq 1 ' ) ] 

< 0J*(r) for 
/lW for r € 1) 


(4.66) 


0; (r) = min{ / 0 i _ oo (r),0 i)+oo (r)} and 0/~(r) = maxjft,-oo(r),0i !+oo (r)}, i 


1,2. (4.67) 


The ’’existence” part of Theorem 14.121 will be covered by the following: 
Theorem 4.16. For all r € (0,1) there holds: 


(z) lim 0i,a = 01.400 (t) and 

a—H-oo 

(*i) lim 0i,« = 01,-00 (t) and 

a—>■—oo 


lim 02 , a = 02,400 

a—>-+oo 

lim 02 ,a = 02,-00 

a—>■—oo 


(4.68) 

(4.69) 


By virtue of (I4.53P and (|2.51l) . we know already that Theorem 14. 161 holds for the Toda system 
r = 1/2; so we only need to establish (14.6811 and (|4.69l) when, 


r € (0,1) and 

We start to analyse (|4.68l) . and then handle (14.69[) similarly. Actually, we are going to 
establish (14.6811 along any sequence: 


where we let (up to a subsequence): 


a. n —> +oo, 


with, 


01 , n ■= 01 («n) ->■ 01 and 02 ,n ■= 02 («n) 02 , 


max{4(A r + 1 ),0 1 (t)} < 0i < 0i and 0 2 = <0/(0i), 

or equivalently, 


max{4,0 2 (r)} < 0 2 < 0 2 and 0i = <//j"(0 2 ). 
Also, to simplify notations, we let: 

vi, n (r) ■= vi(r,a n ) and v 2 , n ■= v 2 (r,a n ), 


(4.70) 


(4.71) 


(4.72) 
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and recall that, 


fi,n(0) = a n -P- +oo, and u 2 , n (0) = 0. 

Therefore, we see that v\ jn admits a blow-up point at the origin, in the sense of Brezis-Merle 
El; see also [ 3 | 68 j for more details. 

Actually, as a first important information, we prove that also v 2) n must admit a blow-up 
point at the origin. 

To this purpose, we recall a well known result established in [3] concerning blow-up sequences 
for ”‘singular’” Liouville equations, which we formulate in a form convenient for our purposes. 
Therefore, for p > 1, we consider v n G W 2,P {B R ) and a n G L P (B R ), satisfying: 

-Av n = \x\ 2N e Vn + a n in B R 

sup R „ v n —> +oo 

Br (4.73) 

(Ib h \x\ 2N e Vn dx) + hn\\ L p( BR ) < C 

sup dBR v n - inf 9Br v n <C 

with N > 0 and a suitable constant C > 0. Assume in addition that v n admits the origin as 
its only blow-up point in B R , in the sense that: 


We have: 


Ve>0 3 > 0 : sup v n < C £ . 

£<\x\<R 


(4.74) 


Theorem 4.17 ([3]). Let N > 0 and p > 1. If v n satisfies (14.731) and (14.741) then, along a 
subsequence, the following holds: 


X \ 2N e Vn —*■ 8n(N + 1)(5q as n —>■ +oo, 


weakly in the sense of measure in B R . 

□ 

Going back to our sequence (i^n, ^ 2 ,n), w e can use Lemma H~T1 and Lemma fOl in order to 
claim the following uniform estimates: 


V r > 0 : vi tn (r) + 2(N + l)logr < C 0 v 2 , n (r) + 2 log r < C 0 , (4.75) 


V 0 < r 0 < -R 0 : sup v-j !n (r) < 7 inf v jin (r) + C(1 + logr 0 ), j = 1 , 2 , (4.76) 

[r 0 ,Jio] fo.-Ro] 

with suitable constants: Cq = Cq(t, N ) > 0, 7 = 7 (-^) G (0,1) and C = C(r, N, > 0. 
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From (14.7511 . we see that Vi >n attains its maximum value at some Ri )Tl € [0, oo), and therefore, 


Vi, n (Ri,n) = max Vi, n (r), i = 1,2, (4.77) 

0<r<oo 

and, by taking also into account (14.7511 . we have: 

vi,n{Ri,n) —> +oo and R\^ n —>• 0, as n —> +oo. (4.78) 

We show that the same property holds for V 2 , n - 

Lemma 4.18. The sequence V 2 n admits a blow up point at the origin, and more precisely: 

V 2 ,n{R 2 ,n) —>• +oo and i? 2 ,n —> 0 , as n —>• +oo. 

Proof. We argue by contradiction and suppose that, for suitable ro > 0 and Cq > 0 we have: 


•y 2 ,n(r) < C 0 , VO < r < r 0 . 

Furthermore, we can take ro > 0 smaller if necessarily, in order to ensure that, 

fro 


(4.79) 


pr 

Jo 


re v 2 ,n(r) < 2. 


To proceed further, we recall the following consequence of the Alexandrov-Bol’s inequality 
(see 0), as pointed out by Suzuki in Proposition 4 of [55] : 

Theorem 4.19. Let 17 C M 2 be a bounded domain with smooth boundary d 17. If p € C 2 (17)n 
C°(17) satisfies 

— Alog p < p in 17, (4.80) 


and H := / p(x)dx < 87 r, then we have: 

Jn 


-2 


( £ Y 

max v < 1 -max p. 

n V 87r / 


(4.81) 


□ 


Clearly, Theorem 14.191 applies for the radial function p = e V2 ’ n in 17 = B ro , and for every 
0 < r < ro, it implies that: 


1 = e V2 ’ n < rnaxe^ 2 ’" < [ 1 — 


B r 


^ fo^ 


-2 


e V2, n (r) < 4 e V2,n (r)_ 


(4.82) 


In other words, by combining (I4.79H and (14.8211 with ([4.7511 and (I4.76H . we conclude that, 


V i? > 0 3 Cr > 0 : \\V2,n\\L°°(B R ) < Cr. 


(4.83) 
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As a consequence of (14.831) . we are in position to apply Theorem 14.171 to the sequence vi iU , 
and along a subsequence, we conclude that: r 2 N e vi ’ n —>• 8 tt(N + l)<5o, weakly in the sense of 
measure on compact subsets of M 2 . 

Furthermore, since e Vl ’ n is uniformly bounded away from the origin (see (14.751) ). by well 
known elliptic estimates, we also deduce that, 


Ve>0 3 C £ > 0 : W v 2 ,n\\c 2 ’ a {B l / e \B e ) < C £ . 

Thus along a subsequence, V 2 , n -T u 2 uniformly in C 2 oc (M 2 \ {0}), with V 2 satisfying, in the 
sense of distributions, the following problem: 


(4.84) 


—Av 2 = e V2 — 8 tt(N + l)ri5o, in Ik 

I ^211 (Bi ) < Ci 

and this is clearly impossible. Therefore, by taking into account also (14.751) . we may conclude 
that: V 2 , n (R 2 ,n) —>• Too and i? 2 , n —>• 0, as n —>• Too, as claimed. 

□ 


Remark 4.20. A crucial ingredient in the proof of Lemma 4 .18, is given by the information 
that V2, n (0) = 0. In fact, we can check easily that the same conclusion would hold under the 
more general assumption: V 2 , n (s n ) = C with s n —>• 0, as n —>• Too. Clearly, the role between 
vi tn and V 2 >n can be interchanged. 

Next, for > 0 defined in (14.771) . we let: 


v l,n (-^1 ,n) 

£1 ,n = e~ 2 ^+D 


v 2,n( R 2,n) 

£2 ,n = e 2 


(4.85) 


so that, from (14.781) and lemma 14.181 we know that, 


We define: 


Ri, n -T 0 and £i >n —> 0, as n —>• 00 ; * = 1,2. 


(4.86) 


v i,n( r ) = v i,n(£2, n r) T 2(N T 1) loge 2 ,n, 

and 

v 2 ,n( r ) = V2,n(£2,nr) T 21og£ 2 ,n- 
We know that, (u* n , n ) satisfies (14.61 p and, 
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r+oo r+oo 

I r 2N + 1 e v i, n ( r ) dr = r 2N+1 e vl ' n ^dr = fi\ n —» /3\ as n —>• oo, 

0 + oo rJ° ( 4 - 87 ) 

as n —>• oo. 


/*~rOO /*H-CO 

/ re V 2,n( r ) dr = re v 2 ,n(r) ^ ^ ^ 

io io 

Furthermore, we notice that: 


maxr>o[wi,„(r)] = v*, n ) = 2(iV + 1) log (g*) , 
maXr>ob2,n( r )] = *>2,n (^) = °- 

We have: 

Lemma 4.21. Along a subsequence, the following holds as n —>• oo : 


-Ri,i 

i? 2 . 


El,?- 


h 1 -!. 


0 , -^-> 0 , 

^2,n ^l,r; 


0 and 


R<2,r 

£ 2 ,n 


(4.88) 

(4.89) 


(4.90) 


with suitable l > 0. In particular v\ n admits a blow up point at the origin while is 
uniformly bounded in C^“(M 2 \ { 0 }). 


Proof. From Lemma [4.181 we derive in particular that R 2 , n > 0, and therefore we can use the 
second equation in (14.6111 to find, 


rRi.n | rr 

0 < V 2 n (R'2,n) ~ V 2 , n (0) = ~ ~ (:te ~ Tt 2N+1 e V ^)dt 

Jo r Jo 


< T 


nr 


r 2(7V+1) 

t 2N+l e vi,n(t) K IIh 2,n 


(4.91) 


(2(N + l )) 2 




Thus, if we recall that: v 2t1l (0) = 0, R 2 ^ +1 ^e v ^-,n(Ri,n) < (j ( see (|4.75l) 4 and v 2:n (R 2tn ) —> +oo 
as n —>• +oo, from (14.9111 we conclude that, 


jf{N+ 1) 

2 (j V+1 ) ~t 0) as n —>• +oo; (4.92) 

^ 2 ,n 

/ t~> \ 2(7V+1) 

( —— ) = ^ 2 ^ + 1 ^_|_qq^ as n + 00 . (4.93) 

V e l ,n ) 

Furthermore, by (14.9311 and the fact that R 2 71 e V2 ’ n( ' Rl ’ n ' > < C (see (14. 7511 ). we also deduce that, 


( £2 ' n 
\£l,n 


2 



^2,n \ 2 _1_ > J_ 

El ,n / (R 2 , n ) 2 e V ^ R ^ ~ C 


U-2,n 
El ,n 


2 


—^ +oo, 


as n —>• +oo. 
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In other words, as n —>• +oo there holds: 


£2 ,n 
£l,n 


+OO, 


and so, 


«l,n f 7 ^) = 2 ( jY + !) log f -> + 00 , 

\ £ 2,n J \ £ l,n J 

< ^^2 ne V2, n (R2,n)/2 < C ^hn ^ 0 

£2,n ^2 ,n ’ -^2,71 

Thus, we have shown that, ujp n admits a blow up point at the origin, and by the estimates 
(14.75p (I4.76|) . which continue to hold for the pair (v$ n ,V 2 n ), we see that actually the origin 
is the only blow-up point for v\ n . 

On the other hand, from (14.8911 we see that v\ n is uniformly bounded from above, and 
v 2 ,n(^ i2 ') = 0- At this point, by recalling that: < C (see (14.751) 1. we can use Remark 

14.201 in order to conclude that must be also bounded from below away from zero. So, 

£2,n 

along a subsequence, there holds: 


—— —» l > 0, as n —>• + 00 . (4.94) 

£2 ,n 

Consequently, we can use (14.9411 together with the Harnack-type inequality (I4.7611 . in order 
to conclude that, for all e > 0 sufficiently small, |u 2 n | is uniformly bounded in [e, ^]. Since 
also v\ n is uniformly bounded from above in [e, j], we can use standard elliptic estimates to 
obtain that v\ n is actually uniformly bounded in Cf ( £( 0 , + 00 ) as claimed. 

Finally, it remains to show that (along a subsequence) tends to zero. To this purpose, 
we recall first that, by (14. 751) . we have: < C. Therefore, if we define: 

Ci ,n( r ) =v* n ( — r) + 2 (N + 1) log f-biir) = ui >n (ei , n r) + 2 (N + 1) log £i,„, 

’ \ £ 2,n J \ £ 2,n J 

£2 ,n(r) = v*2,n (—r) + 2 log f—r) = v 2 , n { £ i , n r) + 2log £ ljU , 

\^2,n J \^2,n J 

we see that (£i,n(r ),£2 ,n(f)) satisfies (14.611) together with the estimates (14.75|) (|4.76|) . 
Furthermore, 


g?2 ,n(r) 

while, £i in satisfies: 



—> 0 uniformly in [ 0 , + 00 ), 
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r > 0 


—A£i ]n = r 2N e^ 1 ' n — re&’ n , 

< max r >o£i,n = (777) = 0 

Jo + °° r 2N+1 e^ 1 ’ n ^ dr < C , 
with a suitable constant C > 0. Therefore, is locally uniformly bounded (see (14.761) 1. and 
by standard elliptic estimates and a diagonalization process, we obtain (along a subsequence): 

1^ £ in C) oc , and - —y Ro? n —y -{-00, (4.95) 

£i ,n 

with £ = £(r) the unique radial solution of the following ’’singular” Liouville equation: 

—A£ = r 2N e £ 

< max r >o6 = ^i(Ro) = 0 
/o + °° r 2N+1 e^ r ' ) dr < C. 

From the classification result in m , we know that, 

£(r) = log ( - 3 (JV+1) ] and f r 2N + 1 f £( r ) dr = 4(1V + 1) (4.96) 

\ 1 + shv+i) 1 ’ ) Jo 

and we can use (14.951) and (14.961) . in order to we derive in particular that Rq = 0, that is: 


as claimed. 


Rl,n 

^1 ,n 


-*■ 0, 


□ 

The proof of Theorem 14.161 

Since n (r) < C, we can use Lemma 14.211 together with Theorem 14.171 in order to conclude 
that, along a subsequence, the following holds : 


r 2N e Vl ’ n —>• 87 t(N + l)5o weakly in the sense of measures on compact sets of R 2 , (4.97) 

and 

v* 2 ,n^v *2 inC 2 c (R 2 \{0}), (4.98) 

as n —>• 00, and defines a (radial) solution (in the sense of distributions) of the following 
problem: 


— Av% = e v z — 87r(JV + 1)tc>o, in R 2 
f R 2 e v Z dx < 00. 


(4.99) 
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By the fact that max r >o u 2 (r) = 0, we can write down the explicit expression for v\ = v%(r) 
in terms of the limit value l in (14.9011 . In particular we know that, 

r+oo 

/ e^ (r) rdr = /3*( r ) = 4 + 8t(N + 1) := $(r), 

Jo 

see [5?j or Theorem 2.2.1 in [6S] ). 

Consequently, from Fatou’s lemma we deduce that, 


r+oo 

lim / e v 2,n( r ) r dr = (3 2 > /3 2 (t). (4.100) 

n^ooj o 

While, by means of Lebesgue dominated convergence, we have: 

Ve>03i? e >0 and n £ € N : 



re 1 ’ 2 ’" 


(■ r) 


dr > P*(t) 




Vi? > i? e and Vn > n e . 


(4.101) 


Claim 1: if r € (0, Tq 1 ^] then the pair (ft, ft) satisfies: ft = 4(iV + 1) and ft = ftft). 

(4.102) 

Indeed for r € (0, Tq 1 ^], we can use (14.4311 to obtain that, 

<pt(PZ) = 4(JV + 1 ) < ft = vtih) < <^2 (ft*), (4.103) 

where, the last inequality follows by (14.1001) and the monotonicity of Clearly, (14.1021) 
easily follows by (14.1031) . 

Claim 2: if r G (tq 1 ^,^ 1 ^) then ft = ft* and ft = ft (4.104) 

To establish (I4.104p . we notice that from (14.481) in this case there holds: 

ft**(r) = p>2 (ft ( r )) > ft“(ft) = ft- (4-105) 

While, from part (i) of Proposition 14.101 we have: 

ft - rft > ft(r) - rft*(r) > 2. (4.106) 

We £x Eo > 0 sufficiently small, so that: 

ft*(T)-T/?r(r)> 2(1+ e 0 ), (4.107) 

and then, we let no £ N and i?o > 1 sufficiently large, so that V n > no there holds, 
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(see (I4.101D ). and 


rfto 

/ re V 2 ,n( r ) dr > ft — £(b 

jo 

roc 

/ r 2iV +Vi*,» dr < ft* + e 0 , 

Jo 


(4.108) 


(4.109) 


(see mm)- 

As a consequence, for every r > Rq and n> n 0 we can estimate: 


v 2,n( r ) - v 2,n( R o) = ~ f R \ (J Q i^ 2 ’^ ~ TS 2N+1 e V l”^) dsj dt 

pr pRo 

= —(logr) / (se V2 ' n ^ — Ts 2N+1 e Vl ' n( ' s ' > ) ds + (log.Ro) / (se V2 ’ n ^ — rs 2N+1 e Vl ’ n ^) ds 

Jo Jo 


+ / (log s)(se V2 ’ n ^ — Ts 2N+1 e Vl ’ n ^) ds. 

Jr 0 


Thus, from (14.1101) . we obtain: 


(4.110) 


rri 0 rv 

ftftr) - v 2 ,n( R o) < -(log -^-)( / se^.nO )ds -t s 2N+1 ds). (4.111) 

Ho Jo Jo 

Therefore, by recalling that v| n (.Ro) < 0, we can use (14.1081) and (14.1091) and (14.1111) to find 
a suitable constant Cq > 0 such that , 


v 2,n( r ) < —(logr , )(/32("T") - rft*{r) - (i + r)e 0 ) + C 0 , for n > n e , r > R 0 , (4.112) 

where, form (14.1071) we know that, 

(ft ( T ) " rft*(r) - (1 + r)e 0 ) > 2 + (1 - r)e 0 > 2. (4.113) 

As a consequence of (14.1121) and (14.1131) . we see that, for every e > 0 we can find R e > 0 and 
n e € N such that, 


p+OO 


re v 2 ,n( r ) dr < —, Vn> 






By taking into account that, 

fRe pRe P+OO 

lim / re v 2 ,rS T 'i dr = re v ' 2 ( ? ) dr < re” 2 ^) dr = ft, 

n_> +oo JO JO Jo 

and (14.1001) . we can take n e larger if necessary, so that Vn > n £ the following holds: 


(4.114) 
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rRe p r+oo 

/ re w 2,n( r ) dr < f3\ + and (3% — £ < / re 112 '"^ dr. 

Jo 2 ' J 0 


In conclusion, for every e > 0, we have shown that, 

r *+00 

/ 0 

and so: 0?j = 02 and 0i = ^(/3|) = (3 as claimed. 


r+oo 4 

02 — £ < / re 1 ’ 2 ’"^ dr < f3 2 + £, Vn> n £ ; 

Jo 


(4.115) 


(4.116) 


Claim 3: if r G [t-^, 1) then 02 = f3 2 and 0i = / 0 1 (4.117) 

The case r = r[ 1 ' 1 follows by a limiting argument form Claim 2. Indeed, asr / t-^, we have: 
02(t ) —t 02 (ti^) = 0 2 (t} 1 )) and 0 **(r) — »• 0i*(t 1 ^) = 0 1 (T-j 1 ^). So it will suffice to prove 
(14.10411 when r G (t['0 1). 

Since in this case: ^(r) — t0**(t) < 2, then the inequality (14.10011 must hold with the strict 
sign, that is: 

02 >0 2 *- (4.118) 

To establish (I4.117H . we need to show that 02, T = 02~ *"01 = 2. To this purpose, we recall that 
02 “ t0i > 2, and so arguing by contradiction, we suppose that actually: 02,r = 02 ~ t0i > 2. 
Let us define (via Kelvin transform): 


Vi,n( r ) = V i,n ( l j + 0t > 10 S f \ ) * = 1,2, 


which satisfies: 


with 


and 


—Avi iU = 2 l Ar+2 ) e l ’ 1 ’ n — rr^,r 4 e i; 2,r 1 

— Av2 ,n = r^ 2 ’ T ~ A e V2 ’ n — r?'4i,T^ 2 ( iV + 2 ) e ^i,n 

0"r = 01,n - 1"02,n and 0£ T = 0 2 , n - T0i,„, 


/*00 7-00 

/ = 0i iW and / s^- 4 e i,2 ^ s Us = 0 2 ,. 

00 ’ 00 


We set, 


and 


0?° := lim(liminf [ 2lyN+2 ' > e vl ’ n ds), 

r —>0 n—>-oo Jq 

rr 

07° := lim(liminf / s^ 2 ’ T ~ A e V2 ' n ds). 

r—^0 n—Kx> Jg 


(4.119) 


(4.120) 


(4.121) 
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Recalling that, V2, n is uniformly bounded in C^“(M 2 \ {0}), while for every compact set 
0 CC M 2 \ {0} we have: maxhi^ —> —oo, we may deduce that, along a subsequence, the 
following holds: 

V 2 ,n —>• V 2 uniformly in C 2 (M 2 \ { 0 }), 

and 

r /h,T~ 2 ( Ar + 2 ) e *h,n —*■ /3^°5q weakly in the sense of measures, locally inK 2 . 
as n —> oo. Furthermore, as /? 2 , r > 2, we have also that, 

r l 3 2 ,T- 4 e v 2 ,n f 3 ™§ 0 _)_ r / 32 ,T- 4 e «2 as n —> oo, 


weakly in the sense of measure on compact set of M 2 . In particular, v 2 defines a radial solution 
for the problem: 


-Av 2 = \x\^~ 4 e i2 + (/3f° - r/3f )<5 0 in M 2 

f R 2 \x\^ T ~ 4 e V2 < 00 . 


(4.122) 


As already mentioned, all solutions of (14.1531) are completely classified in [57] and in partic¬ 
ular, in order to fulfill the integrability condition, we must have: 


/v-(/3 2 °°-T/3n> 2 . 


(4.123) 


To proceed further, we fix £ > 0 sufficiently small and let n £ € N and R e » 1 sufficiently 
large so that, for any n> n £ and R> R e the following holds: 

F F /*°° F 

/ 32 T >/3 2 t - 7 , / se v *^ds < + 7 and / s 2N+1 e v ^ds > ^ (4.124) 

4 Jr 4 J R 4 r 

Therefore, for r > iL and n > n £ we obtain: 


t 


,n( r )~ v 2,n{ R e) = - t [ se v ^ds)dt + t f f s 2N+1 e v ^ds)dt 

Jr, 1 Jo JRr 1 Jo 


ft 

t=r 

rt 

t=r ri 

—logt / se l ' 2 ’ n ds 

+ / flog te V2 ’ n dt T r 

log t / * 2 '' v+1 e ,J * ; n f i. s . 

~T 

Jo 

t=R e J Re 

Jo 

t=R e JR 


< - 


log r (^J se V2 ’ n ds — T J s 2N+1 e Vl ’ n ds\ + J log s (se 1 ’ 2 ’™ — TS 2N+1 e Vl ’ n ^j ds + C £ 

(4.125) 


with C £ > 0 a suitable constant depending on e only. 
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Therefore, by recalling that /3 2T = J 0 °° — rs 2N+1 e Vl ds we find: 


V*2,n( r ) - V2,n(Re) + /? 2 , r log T 

roc 

= logr / se V2 ’ n ds — 


roc roc rr rr 

/ se V2 ’ n ds — rlogr / s 2N+1 e Vl ’ n + / slog se 1 ' 2 ’" — r / s 2Ar+1 log se^^ds + C £ 
Jr Jr J R e J R e 


< log rj se V2 ’ n ds — rlog r / s 2N+1 e Vl ’ n + log r / se V2 ’ n — t 


L 


r 

Jr 


■ f s 2N+1 log 

jRe 

roc roc rr 

= logr / se V2 ’ n ds — rlog r / s 2N+1 e v ^ n — r / g^+ijog S e Vl ’ n ds + C £ 

J Rr Jr j Rr 


se Vl ’ n ds + C e 


< logr 


roc roo 

/ se V2 ’ n ds — t / s 2JV+1 e l ’*> n 
_Ji? e Jr 


+ C'e 


(4.126) 


and by means of (14.1241) . we derive: 

v*2 ,n( 7 ’) < -logr [P 2 ,t ~ (P? ~ r/3 i°) -e\+C £ 

with C e > 0 a suitable constant depending on e only. At this point, we can use (14.1231) . to 
obtain that, for e > 0 sufficiently small, there holds: 


v 2 ,n( r ) < — (logr) (2 + e) + C £ , for r > R e and n>n £ . (4.127) 

But this is impossible, since as above it yields that: f3 2 = P%, in contradiction to (14.1181) . 
Thus we must have that necessarily: p 2 — t/3\ = 2, that is 02 = P 2 ( r ) an d Pi = Pi( T )i and 
also Claim 3 is established. 

Since Claim 1, 2 and 3 hold along any sequence, clearly they imply that part (i) of Theorem 
14.161 holds. 

In order to establish part (ii) we argue similarly, only now we take a sequence: 


a n —>• — oo 

and as above, we let Vi >n (r) = Vi(r,a n ), i = 1,2 and suppose that, 

r+oo p+ oo 

Pl t n := / s ™+l e v 1 :n (s) ds /3] and n / se v 2 ,„(*) ds p 2 

Jo Jo 

with Pi and p 2 satisfying (|4.71[) and (14.721) . 

In view of (14.1281) . it is convenient to let, 

_ a n 

t n = e 2 (iv+i) -poo, as n —>• oo, 

and consider the following rescaled version of (vi tU ,v 2>n ) : 

vi.nW = vi ;n (t n r) - ug n (0) = vi tn (t n r) - a n = vi >n (t n r) + 2(7V + 1) \ogt n , 
V 2 ,n(r) = v 2 , n (t n r) + 21og t n ; 


(4.128) 


(4.129) 
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Thus, (vi >n ,V 2 , n ) is a solution of the system of ODE’s in (14.611) and it satisfies: 


r*+oo 


r»+oo 


S 2AT+I e fi ljn ( s ) ds= S 2N+I e v 1:n (s) dg = _> ( g i 


f‘-\~00 p-\-o o 

/ se V2 ’ n ^ ds= se v 2, "( s ) ds = /3 2in ~+ /?2 

do do 

vi, n (0) = 0 and h 2 , n (0) = 21ogf n ->• +oo 


(4.130) 


as n -> oo. Therefore, the arguments above can be applied to (ui, n , simply with the 
role of v\ )U now played by v 2 ,n and that of v 2>n played by v\, n . 

Hence, with the obvious modifications, one can carry out the same blow-up analysis and as 
above (along a subsequence) arrive at the following conclusion: 


r+o 0 

r+o 0 1 

/ s 27V+1 e Cl ’" (s) ds = 

/ ds -> ^ 1,-00 = ^ 

Jo 

do 1 


Pi(r) 


for r G (0, t| 2 ^) 
for r G [r{ 2 \ 1) 


while, 


r*+oo 


se 


V2,n( s ) ds 



se v2 ’™ (s) ds->¥>+(&) = &,_«> 


4 for r G (0, Tq 2 -*) 

' P 2 *(t) for r G (t,S 2 ) ,t| 2) ) 
P 2 (r) for r G (t 2 (2) ,1) 


as n —> +oo, and this completes the proof of Theorem 14. 161 

□ 

Next, we proceed to show that actually the given condition (14.591) ( or equivalently (14.60(1 ) on 
the pair (/3i,/3 2 ) is also necessary for the solvability of (14.1I) - (I4.2I) . According to (14.541) and 
(I4.55|> . it is clear that we need to be concerned only with the case: r G (0, t-^). 

To this purpose we establish the following: 


Theorem 4.22. 


(*) //tG(0, -) then f3i < (3* i = 1,2, and /3 2 = pt(Pi) > max{/3 2 *,4}, 

or equivalently (3\ = tp^iPz) > max{/3”,4(N + 1)}. 

1 (2) ** * (4-131) 

(ii) If t G (—, ) then fi 2 < /3 2 and Pi > P* > 4 (IV + 1). 

(in) If r G (-jt} 14 ) t/ien /?i < /?** and P 2 > P 2 > 4- 

Proof. It is convenient to introduce the following change of variable: r = e t , and so consider 
the functions: 

z(t) := v\(ef) and u(t) := v 2 (e t ), t G 1; (4.132) 
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satisfying: 


+ e 2(Af+l )t+z _ Te 2t+u _ q for tel 

<pu + e 2t+« _ re 2(iV+i)t+ 2 = Q for f e R 

^ (4.133) 

^j|(—oo) = ^(—oo) = 0, z(—oo) E R, ri(—oo) E R 

4 e 2(N+ ') t+z dt = Pi and f R e 2t+u dt = fa¬ 
in analogy to (13.311) we set: 

/(f) ;= f e 2 ^ N+1)s+z ^ds and g{t) := f e 2s+u ^ds, (4.134) 

J— oo J —OO 

so that, /(t) and g(t) define positive and strictly increasing functions, with /(—oo) = 
g(— oo) = 0 and /(+oo) = fa, g{+ oo) = fa. Furthermore, problem (14.1331) can be writ¬ 
ten equivalently as follows: 


= f(t) - rg(t) for (SR 

~W= 9 (t)-rf(t) for t E R. 


(4.135) 


Moreover, by setting: 

*o(f) == e 2 ^+ 1 )^)+ e 2 ^W-2(lV+l)/(t)-2 5 (t) + i/ 2 (t) + i 5 2 (t)-T/(%(t) Vt E R, 

(4.136) 

i(t) : = e 2(iv+i )t+z(t) _ 2 ( N + !)/(/) + l/ 2 (f) vt € R, (4.137) 

and 

V 2 (t) ■= e 2t+u{t) - 2 g(t) + \g 2 (t) Vt € R, (4.138) 

then the conditions (13.321) and (13.331) can be expressed simply as follows: 

^o(t) = 0, 4'i(t) > 0 and ^(t) > 0 Vt E (—oo,+oo]. (4.139) 

Also notice that, by virtue of (14. 136|) and (14. 139|) . we can derive the following identities: 

0 = e 2 (i\r+i) t +*(t) + e 2 t+u(t) + 1 f{t) _ 4(iV + 1} ) + I 5 (f) ( 5 (f) - 4 ) - t / (t)g(t) 

= e 2 (- ?v + 1 V+ z ( t )-|-i/ ( f ) (/(f)— ( 4 (1V+1)+8 t)) (<7 ( f)—2r/(t)) ( 5 (t)- 4 )+e 2t+u W, Vt E R. 

(4.140) 


At this point, we introduce the function: 


R 0 (t) = 2re 2 ( JV+1 )^W (ff(t)-4) + (/(t)-4(JV+l)) ^^V+^+i/^)^/(t)-(4(iV+l)+8r))^ Vt E R, 

(4.141) 
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and by using (14. 13511 . after straightforward calculations we can check that, 


(t) = -(1 - 2 r ) e 2 (^+ 1 )H-*(t) Q (2t + ^(t) - 4 ) + e 2t+u (A , Vt € M. (4.142) 

Similarly for the function: 

Ri(t) = 2re 2t+ “W (/(t) - (4(1V + 1) + 8r)) + 5 (t) (e 2t+ "W + - 4)( 5 (t) - 2r(4(iV + 1) + 8r))) Vf € M, 

(4.143) 

we have: 

^±(t) = -(1 - 2 T )e 2t+U ^ | e 2(^+i)t+*(t) + I(i + 2r)/(t) (/(t) - (8r + 4(7V + 1)))| . 

(4.144) 

In particular by (14.14211 and (|4.144|1 we obtain: 

- ^(t) = (1 - 2r)eW)^) Q(2r + l)(/(t)(s(t) - 4)) + e 2t + u ^ 

- (1 - 2r)e 2t+uW •|e 2 ( iV+1 ) t+2: ( t ) + 1(1 + 2r)/(t) (/(f) - (8r + 4(7V + 1)))| 

= (1 - 4t 2 ) / e 2(N+i)t+z(t) Q 5 (f)( 5 (f) - 4) + e 2t+u(t) ^ 

_ e 2 t+u(t) ^ e 2(N+i)t+z(t) + !/(*)(/(*) _ (8r + 4(1V + 1)))^ |, Vf € R. (4.145) 

We start to analyse the case where: 0 < r < 1/2, and by contradiction we assume that, 

Pi > Pi = 41V + 1) + 8 r. 

Then, for the function: 

H(t) := e 2 ( Ar +i) t + z ( t ) + !/(*) (/( t ) - (4(iV + 1) + Sr)) = e 2 ^ t+z ^ + - P{) 

(4.146) 

we find that, 

H (Too) > 0 and H' (t) = re 2(7V+1)t+z(t) ( 5 (f) - 4). (4.147) 

Therefore, if we let fi G R the unique value such that: g{t\) = 4, we see that H(t) is decreasing 
for t < ti, and in particular H{t) < 0, for every t <t\. 

While H(t ) is increasing for t > t\, and since H( 00 ) > 0, we find a unique to > t\ where the 
function H vanishes. More precisely, 

if (to) = 0 and in particular: /(to) < P* = 4(iV + 1 ) + 8 r, <7 (to) > 4, (4.148) 
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and moreover, 


H(t) := e 2(tf+i)t+*(t) + I/(t)^/(t) - (4(iV + 1) + 8r)) < 0, Vt < i 0 . (4.149) 

As a consequence of (14.1451) . (14.1391) and (14.1491) we derive: 

to ( 

Ri(to) - Ro(to) = (1 - 4r 2 ) J i e 2 (-^v+ 1 )*+«(t) - 4)) + 

—oo ^ 

- e 2t+u V ^ e 2(^+i)frf*(t) + !/(*)(/(*) - (8r + 4(1V + 1)))^ jdi = 

to 

(1-4r 2 ) J <fe 2 ( Ar+1)t+z ( t) v P 2 (t) - e 2t+ “W//(l) jrif > 0. (4.150) 

— OO 


On the other hand, by (14.1411) . and (14.1481) . we see that, 

i? 0 (t 0 ) := 2re 2(7V+1) * 0+2( * o) ( 5 (i 0 ) - 4^ >0, (4.151) 

while by ([4. 1431) . (I4.140P and ()4. 1481) we drive: 


Ri(t 0 ) < g(to ) (e 2 *°+^°) + i( 5 (t 0 ) - 4)(g(t 0 ) - 2r(4(JV + 1 ) + 8 r))) 

< 9(t 0 ) (e 2 *°+^°) + ^(t 0 ) - 4)( 5 (t 0 ) - 2r/(t 0 ))) = 

- s(io) (e 2(Ar+1)to+ " (to) + \f(t 0 )(f(t 0 ) - (4(7V + 1) + 8 r))) = 0 (4.152) 


and so, 


Ri{to) ~ Ro(t 0 ) < 0 , 


in contradiction with (|4.15U|) . Similarly, for 0<r< 1/2 we can check in an analogous way 
that @2 < P I and ( i ) is established. 


i (2 ) 

Next, we check (ii), so we assume that: r € ( 2 > r i ) and we proceed to show that: /?2 < P 2 * ■ 
To this purpose, as above, we argue by contradiction and suppose that, 

P2 > P?- (4.153) 

Therefore from, (14.501) we have: 

Pi = vl(P2)<<pt(P*2*)=Pl (4-154) 
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and so in this case, H{+ oo) < 0. Thus, for the function H in (I4.146|) now we have that, 

H(t) <0 Vi € M. 


At this point, by recalling that r > from (14.1451) we obtain: 

^i(i) - ^(i) = (1 - 4r 2 ) je 2(7V+1)t+ * W T 2 (f) - e 2t+ “ W F(i)| < 0. (4.155) 

On the other hand, under the given assumption, we see that, 

Ri(+ oo) - i?o(+oo) = |/5 2 (/3 2 - 4)(/3 2 - ft*) - |/3i(/3i - 4(Y + l))(ft - ft) > 0, (4.156) 

and, Ri (—oo) = Ro(— oo) = 0. Therefore in view of (14.1551) we reach again a contradiction, 
and (ii) is established 

Finally, we can check easily that (Hi) follows by similar arguments, and the proof is completed. 

□ 

At this point, Theoren M.121 readily follows as a consequence of Theorem 14.161 which takes 
care about the ’’sufficient” part, while Theorem 14.221 provides the proof of the ’’necessary” 
condition. 
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